A COMPLETE VERSION CAN BE OBTAINED AT http://WWW.ASTRO.PSU.EDU/uSERS/MISAWA/PUB/PaPER/4QHIRES.PS.GZ 

Preprint typeset using LAT^X style emulateapj v. 11/12/01 



SPECTROSCOPIC ANALYSIS OF H I ABSORPTION LINE SYSTEMS IN 40 HIRES QUASARS 1 



Toru Misawa 2 , David Tytler 3,4 , Masanori Iye 5,6 , David Kirkman 3-4 , Nao Suzuki 3-4 , Dan 

lubin 3 ' 4 , and nobunari kashikawa 5 ' 6 

A complete version can be obtained at http://www.astro.psu.edu/users/misawa/pub/Paper/40hires.ps.gz 



o 
o 

a 

>—> 

G\ 

6 



> 

o 
o 
p 

o 
l> 
o 



X 



abstract 

We list and analyze H I absorption lines at redshifts 2 < z < 4 with column density (12 < 
log(Yn i/[cm -2 ]) < 19) in 40 high-resolutional (FWHM = 8.0 km s _1 ) quasar spectra obtained with 
the Keck+HIRES. We de-blend and fit all H 1 lines within 1,000 km s _1 of 86 strong H 1 lines whose 
column densities are log Njji > 15 cm~ 2 . Unlike most prior studies, we use not only Lya but also all 
visible higher Lyman series lines to improve the fitting accuracy. This reveals components near to higher 
column density systems that can not be seen in Lya. We list the Voigt profile fits to the 1339 H 1 com- 
ponents that we found. We examined physical properties of H 1 lines after separating them into several 
sub-samples according to their velocity separation from the quasars, their redshift, column density and 
the S/N ratio of the spectrum. We found two interesting trends for lines with 12 < log(A / n i/[cm -2 ]) 
< 15 which are within 200 - 1000 km s~ 4 of systems with log(An i/[cm~ 2 ]) > 15. First, their column 
density distribution becomes steeper, meaning relatively fewer high column density lines, at z < 2.9. 
Second, their column density distribution also becomes steeper and their line width becomes broader by 
about 2-3 km s _1 when they are within 5,000 km s _1 of their quasar. 

Subject headings: quasars: absorption lines — quasars: general — intergalactic medium 



1. INTRODUCTION 

Quasar absorption systems have historically been di- 
vided largely into three physically distinct categories: (i) 
absorption systems with strong metal lines arising in or 
near intervening galaxies, (ii) weak H 1 systems in the Lya 
forest that come from the intergalactic medium, and (iii) 
intrinsic systems that are physically related to the quasars, 
including associated and broad absorption line systems. 

Metal absorption systems usually contain H 1 lines 
with relatively large column densities, including two sub- 
categories: damped Lyman alpha (DLA) system and Ly- 
man limit system (LLS) with H 1 column densities of 
log(A H i/[cm~ 2 ]) > 20.2 and 17.16, respectively. Deep 
imaging observations around quasars have provided ev- 
idences that the metal absorption systems are often pro- 
duced in intervening galaxies. Galaxies have been detected 
that can explain Mg 11 absorption lines (e.g., Bergeron 
& Boisse 1991), and C iv absorption lines (e.g., Chen, 
Lanzetta, & Webb 2001). 

On the other hand, nearly all H I lines have smaller col- 
umn densities (log Nhi < 15) than those associated found 
in gas that shows strong metal lines. The number of H I 
lines per unit z increases with redshift (Peterson 1978; 
Weymann et al. 1998a; Bechtold 1994), because the in- 
tergalactic medium is denser and less ionized at higher 
z. The Lya absorption lines are produced in intergalac- 
tic clouds (e.g., Sargent et al. 1980; Melott 1980), which 
are the higher density regions in the inter-galactic medium 



(IGM). The Lya lines are broadened by Hubble flow (e.g., 
Rauch 1998; Kim et al. 2002a) as well as the Doppler 
broadening from the gas temperature. 

Misawa et al. (2004) presented a study of H I absorp- 
tion lines seen towards 40 quasars in spectra from the Keck 
HIRES spectrograph. In a departure from prior work, they 
considered the H I lines without considering the metal 
lines. They classified the H I lines as either high den- 
sity lines [HDLs] which have or are near to strong H I lines 
that are probably related to galaxies, and low density lines 
[LDLs] that are far from any strong H I lines and are more 
likely to be far from galaxies and hence in the IGM. 

Following Misawa et al. (2004), we define HDLs as all 
H I lines within ± 200 km s _1 of a lines with 15 < log Njji 

< 19 cm~ 2 . We define LDLs as lines with 12 < log Nhi 

< 15 which are within 200 - 1000 km s" 1 of a line with 
15 < log Nhi < 19 cm -2 . This last velocity constraint 
is intended to make the LDL a "control" sample for the 
HDLs, where the two samples come from similar redshifts 
and regions of the spectra with similar signal to noise. 

Misawa et al. (2004) discovered that the HDLs have 
smaller Doppler parameters (b- values) , for a given column 
density than the LDLs, and they also found the same effect 
in a hydrodynamic simulation with 2.7 kpc cells. Misawa 
et al. (2004) suggested that the LDLs are cool or shock- 
heated diffuse intergalactic gas, and that the HDLs are 
cooler dense gas near to galaxies. 

Misawa et al. (2004) fit all the accessible transitions in 
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the H I Lyman series to help de-blend H I lines. Their 
main sample comprised 86 H I absorption systems each 
with log Njji > 15 cm~ 2 . They also fit all H I lines within 
±1,000 km s _1 of these H I lines, to give a total sam- 
ple of 1339 H I lines, including the 86 lines. This is the 
only large sample in which multiple Lyman series lines are 
fit together, although this method has been used on indi- 
vidual systems and small samples (Songaila et al. 1994; 
Tytler, Fan, & Buries 1996; Wampler et al. 1996; Car- 
swell et al. 1996; Songaila et al. 1997; Buries & Tytler 
1998a,1998b, Buries, Kirkman, & Tytler 1999; Kirkman 
et al. 2000; O'Meara et al. 2001; Kirkman et al. 2003; 
Kim et al. 2002b; Janknecht et al. 2006). 

In this paper, we present measurements of the absorp- 
tion lines that Misawa et al. (2004) have analyzed. We 
give a detailed description of each absorption system and 
we summarize new results. The paper is organized as fol- 
lows: In we give descriptions of the data and the line 
fitting. The results of our statistical analyses are presented 
in |J3] We discuss our results in S}4] and summarize them 
in iJ5J In the Appendix, we describe the properties and we 
give velocity plots for each H I system. We use a cosmol- 
ogy in which Ho = 72 km s _1 Mpc _1 , fi TO = 0.3, and f^A 
= 0.7. 

2. SPECTRA AND LINE FITTING 

The 40 quasars in our sample have either DLA systems 
or LLSs, and they were observed as a part of a survey for 
measurements of the deuterium to hydrogen (D/H) abun- 
dance ratio. The detailed description of the absorption 
and data reduction are presented in Misawa et al. (2004) . 
We caution that our sample was biased in subtle ways by 
the selection of LLS and DLAs that seemed more likely to 
show D, i.e. those with simpler velocity structure. 

We list the 40 quasars in Table 1. Column (1) is the 
quasar name, (2) the emission redshift. Columns (3) and 
(4) are the optical magnitude in the V and R bands. The 
lower and upper wavelength limits of the spectra are pre- 
sented in columns (5) and (6). Column (7) gives the S/N 
ratio at the center of the spectrum. Same data set was 
also used in Misawa et al. (2007) in a study of the quasar 
intrinsic absorption lines. 

We will discuss only H I lines with logNni > 15 cm~ 2 
and other H I weaker lines within ±1,000 km s _1 of these 
strong H I lines. We selected this velocity range since it 
is enough to include the conspicuous clustering of strong 
metal lines. Indeed such strong metal lines are normally 
confined to an interval of < 400 km s _1 even for DLA 
systems (Lu et al. 1996b). 

Here we briefly review the line detection and fitting pro- 
cedures that we discussed in more detail in Misawa et 
al. (2004). We began searching the literature for H I 
lines with log N hi > 15, including the DLA and LLS cat- 
alogues (Sargent, Steidel, & Boksenberg 1989, hereafter 
SSB; Lanzetta 1991; Tytler 1982), and metal absorption 
systems (Peroux et al. 2001; Storrie-Lombardi et al. 1996; 
Petitjean, Rauch, & Carswell 1994; Lu et al. 1993; Steidel 
& Sargent 1992; Lanzetta et al. 1991; Barthel, Tytler, & 
Thomson 1990; Steidel 1990a, b; Sargent, Boksenberg, & 
Steidel 1988, hereafter SBS; SSB). We also search for them 
ourselves. If more than one strong H I line was detected in 
a single 2000 km s _1 velocity window, we take the position 



of the H I line with the largest column density (hereafter 
the "main component") as system center. We found 86 H I 
systems with log Nhi > 15, at 2.1 < z a f, s < 4.0, in 31 of 
the 40 quasars. Figure 1 of Misawa et al. (2004) gives the 
velocity plot of one of these systems, and below we give 
the rest. 

We give parameters describing these 86 systems in Ta- 
ble 2. In successive columns list (1) the name of the quasar; 
(2) the redshift of the main component, that with the 
largest column density ; (3) the H I column density of 
the main component Ni\ (4) N2, the second largest H 1 
column density within ±1000 km s _1 of the main compo- 
nent; (5) the ratio of N 2 to Ni; (6) - (9) the S/N ratios 
at Lya, Ly/3, Lye, and LylO; (10) the number of lines in 
the ± 1,000 km s _1 window; (11) the number of H 1 lines 
classified as HDLs (described later); (12) comments on the 
H I system; (13) references. We will call this list sample 
SO (Table 3). 

When we were fitting the lines, we rejected narrow lines 
with Doppler parameter of b < 4.8 km s , which corre- 
sponds to the resolution of our spectra. We also identify 
all lines with b < 15 km s _1 as possible metal lines (called 
M 1 in Tables) because H 1 lines with this narrow width are 
rare (e.g., Hu et al. 1995, hereafter H95; Lu et al. 1996a, 
hereafter L96; Kirkman & Tytler 1997a, hereafter KT97). 
If there was more than one way to fit the lines, we chose the 
fit with the fewest lines. If the model did not give good 
fits to all the Lyman series lines, we adopted the model 
that best fit the lower order lines where the SNR is best. 
For H I lines with column densities of log Nhi > 16.6 the 
Lyman continuum optical depth is r > 0.25. For these sys- 
tems we checked if the residual flux at the Lyman limit was 
consistent. Our fitting method could readily overestimate 
the Doppler parameter but not the column density. Once 
the fitting model is chosen, we used \ 2 minimization in a 
code written by David Kirkman, to get the best fit param- 
eters of H 1 column density (log Nhi), Doppler parameter 
(b), and absorption redshift (z). The internal errors are 
typically c(log N H i)=0-09 cm" 2 , a(b)=2.l km s -1 , and 
cr(z)=2.5xl0- 5 . 

We prepared a sample SI that is a sub-sample of SO in- 
cluding only 973 H 1 lines and 61 H 1 systems with log Nhi 

> 15. SI excludes 25 systems with difficulties such as (i) 
poor fitting due to gaps in the echelle formatted spectra, 
(ii) poor fitting due to strong DLA wings (i.e., log Nhi 

> 19), (iii) close proximity in redshift to the background 
quasars (i.e., within 1,000 km s _1 of the emission redshift), 
and (iv) overlapping with other H 1 systems. The S/N ra- 
tios of the spectra are at least S/N ~ 11 per 2.1 km s _1 
pixel and the mean value is S/N ~ 47 for Lya lines. Among 
these 61 H 1 systems, three systems may be physically as- 
sociated with the quasars based on the partial coverage 
analysis for the corresponding metal absorption lines (Mi- 
sawa et al. 2007). However, we keep these systems in SI 
sample, because we still cannot reject the idea that they 
are intervening systems. 

We give detailed descriptions of all the lines that we fit 
in the Appendix. We also give velocity plots of the first 
five Lyman transitions, Lya, Ly/3, Ly7, Ly<5, and Lye. 
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We investigate the properties of line parameters such as 
the column density, Doppler parameter, and the clustering 
properties of the H I lines. Since this sample contains not 
only H I lines originating in the intergalactic diffuse gas 
clouds (i.e., LDLs), but also H I lines produced by inter- 
vening galaxies (i.e., HDLs), we also attempted to sepa- 
rate H I lines into HDLs and LDLs based on the clustering 
trend (Misawa et a. 2004). 

Our analysis is similar to that of previous studies (e.g., 
H95; L96; KT97), but with three key differences: (i) earlier 
studies used all H I lines detected in the quasar spectra, 
whereas we use only H I lines within ± 1,000 km s _1 of 
the main components with log Nhi > 15, (ii) our sample 
contains a number of strong lines (log Nhi > 15) in addi- 
tion to weak H I lines (log Nhi < 15), and (hi) our sample 
covers a wide redshift range: 2.0 < z a b s < 4.0. The red- 
shift distributions of the 86 and 61 H I systems in samples 
SO and SI are shown in Figure 1. 

3.1. Sub-Samples for the Statistical Analysis 

For further investigation, we prepared several sub- 
samples as follow. It is known that the comoving num- 
ber densities of low-ionization lines, such as H I lines, de- 
creases in the vicinity of quasars (Carswell et al. 1982; 
Murdoch et al. 1986; Tytler 1987). This trend is known 
as the "proximity effect", and is probably caused by the 
enhanced UV flux from the quasar towards which the ab- 
sorption is seen. We separate the 61 H I systems (sample 
SI) into sub-samples S2a (the velocity difference from the 
quasar, Av > 5,000 km s" 1 ) and S2b (Av < 5,000 km s" 1 ). 
We have already removed from SI all H I systems within 
1,000 km s _1 of the quasars to avoid H I systems from the 
quasar host galaxies. 

H95 emphasized that the line detection limit is almost 
wholly determined by the line blending (or blanketing), 
and not by the S/N ratio of the spectrum. In order to con- 
firm whether the distribution of line parameters is affected 
by the quality of the spectrum, we made two overlapping 
samples from SI using the S/N ratio of each spectrum in 
the Lya region: S3a (S/N > 40), and S3b (S/N > 70). 
These sub-samples include 34 (~ 60%) and 17 (~ 30%) of 
the 61 H I systems of the sample SI. 

Sample SI covers a broader range of redshifts, 2.0 < z < 
4.0, when compared with previous studies: 2.55 < z < 
3.19 for H95, 3.43 < z < 4.13 for L96, and 2.43 < z < 3.05 
for KT97. When we investigate the redshift evolution of 
H I absorbers, we also divided SI into two sub-samples; S4a 
(z < 2.9) and S4b (z > 2.9). Here the two sub-samples 
have nearly the same number of H I systems. 

Finally, we made sub-samples according to the column 
densities of H I lines, as the distributional trends of strong 
and weak H I lines are very different (see Figure 2 in Mi- 
sawa et al. 2004); the H I lines with relatively large col- 
umn densities tend to cluster around the main compo- 
nents, while the number of weak H I lines decreases near 
the center of H I systems because of line blanketing. Since 
one of our interests is to determine the boundary value of 
column density between HDLs and LDLs (although other 
parameters may be necessary to separate them), we sepa- 
rate the 973 H I lines into eight sub-samples according to 
their column densities. We use boundary values of log N 
= 13, 14, 15, and 16, where sub-sample S5 a t contains H I 



lines whose column densities are log(A f H i/[cm" 2 ]) values 
of a to b. 

In Table 3 we summarize these 16 sub-samples. Here 
we should emphasize that S5i2i3, S5i2i4, S5i2i5, S5i2i6, 
S5i3i9, S5i4ig, S5i5ig, and S5i6ig, are separated by the 
properties of lines, while S2a, S2b, S3a, S3b, S4a, and S4b 
are separated by the unit of system. 

3.2. Physical Properties o/H I Absorbers 

For each sub-sample prepared in the last section, we per- 
form statistical analysis, including analysis of the column 
density distribution, Doppler parameter distribution, and 
line clustering properties. The samples used here contain 
both HDLs and LDLs. 

3.2.1. Column Density Distribution 

The column density distribution of H I lines, dn(N)/dN, 
are usually fit with a power law (Carswell et al. 1984, 
Tytler 1987, Petitjean et al. 1993; H95), 
, fdn(N)\ 

g (dN ) = ~ 13 gN+ ' (1) 
where the index (3 was estimated to be 1.46 (H95), 1.55 
(L96), and 1.5 (KT97) with only weaker H I lines with 
\ogN H i = 12 - 14.5. Janknecht et al. (2006) found (3 
= 1.60±0.03 at lower redshift z < 1.9. The column den- 
sity distributions of our seven sub-samples (SI, S2a, S2b, 
S3a, S3b, S4a and S4b) per unit redshift and unit column 
density are analyzed. The plot in Figure 2 is the result 
for sub-sample SI. Since the turn-over of the distribution 
around log Nhi = 12.5 is probably due to line blending 
and/or blanketing as described later in § 4, we fit the col- 
umn density distributions only for log Nhi > 13. The 
best-fitting parameters, (3 and A, as well as the redshift 
bandpass, Az, of each sub-sample are summarized in Ta- 
ble 4, along with the past results from KT97 and Petitjean 
et al. (1993). 

The indices that we find, (3 — 1.40±0.03, are slightly 
smaller than the value in the past results, (3 = 1.46 -1.55, 
which means that our sample favors strong H I lines. But 
we expect this type of trend. Our samples contain not only 
LDLs but also HDLs, and since we cover only the velocity 
regions within ± 1,000 km s _1 of the main components, 
we have a strong excess of strong H I lines. We see the col- 
umn density distribution does not change with the velocity 
distance from the quasars (S2a and S2b) or with the S/N 
ratio (S3a and S3b), but it is weakly affected by redshift 
(S4a and S4b). We also applied the Kolmogorov-Smirnov 
(K-S) test to the sub-samples. The results in Table 5 show 
that we can not rule out the hypothesis that they are ran- 
dom samplings from the same population. 

3.2.2. Doppler Parameter Distribution 

The distribution of the Doppler parameter of H I lines 
have been approximately given by the truncated Gaussian 
distribution (H95; L96), 



dn(b) 
db 




b ^ bmin 
b <C bmin 



(2) 



where bo and crj, are the mean and the dispersion of b dis- 
tribution and b m in is the minimum b value for an H I line. 
There are two origins of line broadening: thermal broad- 
ening (&t), and micro-turbulence (btur)- The total amount 
of broadening is given by b = \jb\ + bf ur . 
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In order to determine the correct Doppler parameter, we 
have to individually resolve and fit H I lines using Voigt 
profiles. However, most of weak H I lines disappear in the 
observed spectrum due to line blending and blanketing, 
especially near strong lines. To derive the intrinsic (as 
opposed to observed) distribution of the Doppler param- 
eter, previous authors (i.e., H95; L96; KT97) have cho- 
sen artificial H I lines (as input data) with distributions 
characterized by a Gaussian, and used them for compari- 
son with the observed distributions. KT97 noted that the 
b- value distribution of lines seen in simulated spectra re- 
sembles the distribution of the input data, except for two 
differences: (i) an excess of lines with large b- values is seen 
in the recovered data, compared to the input data, which 
is probably produced by line blending, and (ii) lines with 
very small b- values (b < 15 km s -1 ) are found, which are 
probably data defects or noise, as they are not present in 
the input data. Nonetheless, the 6-value distribution of 
the input and recovered data resemble each other between 
b = 20 km s _1 and b — 60 km s _1 . 

We would ideally like to determine the real distribu- 
tion of Doppler parameters; however, the only way to do 
this is to perform simulations, and compare the recovered 
Doppler parameter distribution with the observed distri- 
bution. Such simulations are expensive to perform, so 
in this work we simply compared the distribution in our 
sample with the past results of H95 and L96 at b = 20 
- 60 km s _1 . We analyzed the observed distributions of 
Doppler parameters for 15 sub-samples, and compare them 
with the results in H95 (bo — 28 km s _1 , = 10 km s _1 , 
and b m in = 20 km s -1 ) and L96 (bo = 23 km s _1 , o-f, = 
8 km s , and b m i n = 15 km s _1 ), where the parameters 
are the inputs used for artificial spectra that reproduce the 
observed distribution. 

We see an excess of H I lines with large 6-values of b 

> 50 km s _1 in all of the sub-samples, while we have no 
lines with b < 15 km s _1 because we decided to classify 
them into metal lines. All sub-samples except for S1213 
have relatively large 6-values, and their distributions are 
closer to H95's distribution than L96's. In contrast, the 
distribution of S'1213, containing only H I lines with small 
column densities, log Nhi < 13, resembles L96's distribu- 
tion. We plot in Figure 3 the Doppler parameter distri- 
butions of sub-samples SI and S1213. We also applied a 
K-S test to the sub-samples (S2a, S2b, S3a, S3b, S4a, and 
S4b). The results are listed in Table 6. The probability, 
that the distributions of sub-samples S2a (Av(z em — z a b s ) 

> 5000 km s" 1 ) and S2b (Av(z em ~ z abs ) < 5000 km s" 1 ) 
were drawn from the same parent population, is very small, 
^2.4 %. This result could suggest that the Doppler pa- 
rameters of H I lines within 5000 km s _1 of quasars are 
affected by UV flux from the quasars. 

3.3. HDLs and LDLs 

In the previous section, we carried out statistical anal- 
ysis using sub-samples containing both HDLs and LDLs 
together. Here, we repeat these tests on the two samples 
separately 

Metal absorption lines seen in DLA systems or LLSs are 
strongly clustered within several hundred km s _1 , which 
implies their relationship to galaxies. In simulations Dave 
et al. (1999) also noted that galaxies tend to lie near 



the dense regions that are responsible for strong H 1 lines. 
On the other hand, for weak H 1 lines, no strong cluster- 
ing is seen (e.g., Rauch et al. 1992; L96; KT97), although 
some studies found only weak clustering trends (e.g., Webb 
1987; H95; Cristiani et al. 1997). 

As Misawa et al. (2004) found, lines with log N HI > 15 
show strong clustering trends at Av < 200 km s _1 , while 
lines with lower column densities cluster weakly at Av < 
100 km s _1 (Figure 4). 

Misawa et al. (2004) defined HDLs as H 1 lines with 
15 < log Nhi < 19 and other weaker H 1 lines within ± 
200 km s _1 of those stronger H 1 lines. They then defined 
the LDLs as all other lines with 12 < log Nhi < 15. They 
chose log Nhi = 15 cm~ 2 in the definition because the 
two point correlation was largest for a sub-sample of H 1 
lines with 15 < log Nhi < 19. We list the number of HDLs 
and LDLs in the sub-samples in Table 7. 

3.3.1. Column Density Distribution 

In Table 8 we give the parameters that describe the 
column density distributions of HDLs and LDLs for five 
sub-samples (SI, S2a, S2b, S4a, and S4b). The most obvi- 
ous result is that the HDLs have a smaller index than the 
LDLs. We see the same result in Figure 3 of Petit jean et 
al. (1993), and hence we now confirm this with the first 
large sample to consider the sub-components of the HDLs. 

The distributions of LDLs in sample SI, S2a, and S4b 
are almost consistent with the previous result in KT97. On 
the other hand, the power law indices for the LDLs of S2b 
(Av < 5000 km s" 1 ) and S4a (z < 2.9), (3 = 1.90±0.16 
and 1.71±0.06, are larger than the values for the other 
sub-samples, (3 ~ 1.52, which means that LDLs at lower 
rcdshift or near the quasars tend to have lower column den- 
sities compared with those at higher rcdshift or far from 
the quasars. The change in the column density distribu- 
tion near to the quasars may be just a consequence of the 
enhanced UV radiation. 

3.3.2. Doppler Parameter Distribution 

The Doppler parameter distributions of HDLs and LDLs 
for sub-samples (SI, S2a, S2b, S4a, and S4b) are also in- 
vestigated, and the results of K-S tests applied to them 
are summarized in Table 9. The only remarkable result 
is that the probability, that the Doppler parameter distri- 
butions of LDLs at Av > 5,000 km s" 1 (S2a) and Ai; < 
5,000 km s _1 (S2b) from the quasars were drawn from the 
same parent population, is very small, ~2.0%. We shows 
these two distributions in Figure 5. In Figure 6 we see 
that the cumulative distribution for the line 6-values rises 
more slowly for the LDLs near to the quasars (at Av < 
5000 km s _1 ), which means that these lines near to the 
quasars are broader by about 2-3 km s . For other pairs 
of the sub-samples, we could not rule out the hypothesis 
that their parent populations are same. 

4. DISCUSSION 

In this section, we discuss our results, especially the fact 
that the column density distribution is changing with the 
rcdshift and the velocity distance from the quasars. We 
also compare our results to those at lower redshift (z < 
0.4) from the literature. After that, we also briefly discuss 
the completeness of H I lines in our 40 spectra. 
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4.1. Redshift Evolution of H I Absorbers 

In § 3, we prepared two sub-samples, S4a and S4b, to 
compare the physical properties of H I lines at lower red- 
shift (z a bs < 2.9) and at higher redshift (z a bs > 2.9). We 
do not see a change in the column density distribution in 
the sample as a whole, but once they are separated into 
HDLs and LDLs, we notice that the index of the column 
density distribution of LDLs at z abs < 2.9 ({3 = 1.71±0.06) 
is clearly different from that of LDLs at z a b s > 2.9 ((3 = 
1.52±0.09). On the other hand, there was no redshift evo- 
lution for HDLs. This trend, shown in Figure 7, means 
that there is a deficit of relatively stronger LDLs (i.e., 
log Nhi > 14.5) at lower redshift. One of the possible 
explanations is that at lower redshift, more H I lines with 
the column densities just below log Nhi — 15 (i.e., rela- 
tively strong LDLs) might be associated with HDLs. In 
other words, stronger (i.e., log Nhi — 14.5 - 15) LDLs get 
into within 200 km s _1 of the nearest HDLs, and would 
be classified into HDLs, which is consistent with the trend 
expected in the hierarchical clustering model (Figure 8). 

As for the Doppler parameter distribution, we did not 
find any remarkable redshift evolutions in neither HDLs 
nor LDLs. L96 claimed that there is a redshift evolu- 
tion of the Doppler parameter between z = 2.8 and 3.7; 
the mean value of Doppler parameter at z a b s = 3.7 (6o 
= 23 km s _1 ; L96) is smaller than the value at z a bs = 
2.8 (bo — 28 km s _1 ; H95). The corresponding value in 
KT97 (bo = 23 km s _1 ) is, however, different from the 
result in H95. The difference may be due to the differ- 
ent line fitting procedure used in these studies; L96 and 
KT97 used the VPFIT software, while H95 used different 
software. Especially important is how the authors chose 
to treat blended lines. The difference could be related to 
the difference of the spectrum resolutions; R=45000 (L96; 
KT97) and R=36000 (H95). Janknecht et al. (2006) did 
not detect any evolution on the Doppler parameter at z 
= 0.5 - 1.9. Our results, which are based on the data set 
taken with one observational configuration and fit using 
the same procedure, suggests that the Doppler parameter 
distribution of H I clouds does not evolve with redshift at 
z = 2-4. 

4.2. Proximity Effect near Quasars 

It has long been noted that the number of Lya lines 
decreases near to the redshift of the quasars (Carswell et 
al. 1982; Murdoch et al. 1986; Tytler 1987). This phe- 
nomenon is related to the local excess of UV flux from 
the quasars. The proximity effect has been used to evalu- 
ate the intensity of the background UV flux. Bajtlik et al. 
(1988) first measured the mean intensity of the background 
UV intensity, J v = l()- 21 - 0±0 - 5 (erg s _1 cm" 2 Hz" 1 str" 1 ) 
at the Lyman limit at 1.7 < z < 3.8, by estimating the 
distance from the quasar at which the quasar flux is equal 
to the background UV flux. The typical radius is ~ 5 Mpc 
in physical scale that corresponds to the velocity shift of 
Av <~ 4,000 km s _1 from the quasars. L96 also evaluated 
the background UV intensity to be J u = 2 x 10~ 22 (erg s _1 
cm" 2 Hz" 1 str" 1 ) at z ~ 4.1 in the spectrum of Q0000-26. 

We see two differences in LDLs within 5,000 km s _1 
of the quasars, compared to those far from the quasars 
(Av > 5,000 km s^ 1 ). We see fewer strong LDLs leading 
to a large index for the column density power law, (3 = 



1.90±0.16 (Figure 9). We also see that the distribution of 
Doppler parameter is different from that of H I lines far 
from the quasars at 98.8 % confidence level. The lines near 
to the quasar apparently tend to have broader lines (Fig- 
ures 5 and 6), although this is a tentative result because 
we consider few lines near to the quasars. 

These results could be accounted for by assuming a two- 
phase structure: outer cold low-density regions and inner 
hot high-density regions in which temperature is deter- 
mined by the competition between photoionization heat- 
ing and adiabatic cooling. When gas is near to the quasars, 
the outer regions become too highly ionized to show much 
H i, and only the inner hot regions would be observed in 
H i, which would increase the mean value of the Doppler 
parameter. The increased ionization also decreases the 
total column densities of H i gas compared with gas far 
from the quasars (Figure 10). As reported in the past ob- 
servations (e.g., Kim et al. 2001; Misawa et al. 2004), 
log Nhi and 6(H i) have a positive correlation for log Nhi 
< 15. This correlations was also reproduced by hydro- 
dynamical simulations (e.g., Zhang et al. 1997; Misawa 
et al. 2004). These results suggest that high density re- 
gions tend to have larger Doppler parameters, if the ab- 
sorbers are not optically thick. Dave et al. (1999) also 
presented an interesting plot in their Figure 11 that sup- 
ported there existed three kinds of phases for H I absorbers 
(diffuse, shocked, and condensed phases). Among them, 
the diffuse phase whose volume densities are small (i.e., 
it corresponds to LDLs in our paper) has a positive cor- 
relation between log Nhi and b. On the other hand, an 
anti-correlation between log Nhi and b is seen only for the 
condensed phase with high volume density that is proba- 
bly associated with galaxies. The shocked phase, probably 
consisting of shock-heated gas in galaxies, does not show 
any remarkable correlations between them. Thus, if we 
assume all LDLs in our sample arise in the diffuse phase 
absorbers, our scenario above could reproduce the differ- 
ence between sub-samples S2a and S2b. 

4.3. Comparison to H I Absorbers at Lower Redshift 

The number density evolution of H I absorbers (i.e., 
dN/dz oc (1 + z)~ 7 ) has been known to slow dramati- 
cally at z ~ 1.6, from a high-z rapid evolution with 7 of 
1.85±0.27 (Bechtold 1994) to a low- z slow evolution with 
7 of 0.16±0.16 (Weymann et al. 1998b). This trend is sug- 
gested to be due to the decline in the cxtragalactic back- 
ground radiation using hydrodynamic cosmological simu- 
lations (e.g., Theuns, Leonard, & Efstathiou 1998). Thus, 
a comparison of H 1 absorbers at high-z and local universe 
is another interesting topic. 

In § 3, we found that the column density distribution of 
LDLs at z < 2.9 (f3 = 1.71±0.06) is steeper than that at 
z > 2.9 ((3 = 1.52±0.09). We proposed this trend could 
be due to the hierarchical clustering. If the assembly of 
structure in the IGM indeed dominates the column den- 
sity distribution, we would expect to find a steeper column 
density distributions at lower redshift as proposed in § 4.1. 

Using the Hubble Space Telescope (HST) and the Far Ul- 
traviolet Spectroscopic Explorer (FUSE), Penton, Stocke, 
& Shull (2004) and Lehner et al. (2007) estimated power- 
law indices to be (3 of 1.65±0.07 at 12.3 < logN HI < 14.5 
and 1.76±0.06 at 13.2 < log N HI < 16.5 at z < 0.4, respec- 
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tively. On the other hand, Dave & Tripp (2001) found a 
flatter distribution (J3 = 2.04±0.23) at z < 0.3. The latter 
steeper distribution was also reproduced by hydrodynami- 
cal simulations (e.g., Theuns et al. 1998). If we accept the 
steeper result, the column density distribution would con- 
tinue to be steeper as going to the lower redshift, which 
supports our idea that the hierarchical clustering could 
play a main role of the evolution seen in the column den- 
sity distribution, although extragalactic radiation would 
contribute to play a role. 

Absorption line width is another parameter that is still 
in argument whether it would evolve with redshift or not, 
as mentioned in § 4.1. While most of the space-based ul- 
traviolet observations could not measure line widths by 
model fittings because of the lacks of spectral resolutions, 
Lchner et al. (2007) for the first time measured Doppler 
parameters of H I absorption lines accurately at lower red- 
shift (z < 0.4), and investigate their distributional trend. 
By comparing to the results at higher redshift, Lehncr 
et al. (2007) discovered that Doppler parameters are 
monotonously increasing from z = 3.1 to ~0. Such a 
trend was not confirmed in past papers (e.g., Janknecht 
et al. 2006). The fraction of the broad Lya absorbers 
(BLA; b > 40 km s _1 ) is also confirmed to increase by a 
factor of ~3 from z ~ 3 to (Lehner et al. 2007). Here, 
b = 40 km s _1 corresponds to gas temperature of T gas ~ 
10 5 K, which is a border between the cool photoionized ab- 
sorbers and the highly ionized warm-hot absorbers. These 
results suggests that a large fraction of H I absorbers at 
very low redshift (i.e., z < 0.4) are hotter and/or more 
kincmatically disturbed than at higher redshift (i.e., z > 
2.0). 

In our sample, we do not see any clear difference of 
mean/median Doppler parameter at z > 2.9 (b mean — 
31.0±10.0, b med = 28.1) and z < 2.9 (b mean = 32.0±10.9, 
bmed = 29.6). Neither HDLs nor LDLs shows any evo- 
lutional trends. These negative results could be because 
with our optical data we covered only higher redshift re- 
gions than z ~ 1.6, at which dN/dz evolution dramatically 
changed. Similarly, the fraction of the BLA (Jbla — 0.182 
at z > 2.9 and 0.196 at z < 2.9) shows an only marginal 
hint to the evolution. However, these fractions are con- 
sistent to the result from KT97 (Jbla — 0.179; Lehncr 
et al. 2007) at 2.43 < z < 3.05 that is similar redshift 
coverage as our sample. Thus, Doppler parameter could 
increase as going to the lower redshift, but such a trend 
would be remarkable only if we trace its distribution at 
very low redshift (at z < 0.4) and compare it to that at 
much higher redshift (at z > 2). 

As for the clustering trend of H I absorption lines, we see 
a very similar property at low and high redshift regions. 
As presented in Figure 4, we found a strong clustering 
trend within An of 200 km s _1 for H I lines with log Njji 
between 15 and 19, while only a weak correlation is seen 
for weaker H I lines within Av of 100 km s _1 . Penton et 
al. (2004) presented very similar results: 5ct (7.2c) excess 
within Av of 190 km s" 1 (260 km s" 1 ) and only stronger 
H I lines contribute to this clustering. Penton, Stock, & 
Shull (2002) proposed such clustering trends within sev- 
eral hundreds of km s _1 are due to clusters of galaxies. 
There could exist similar kinematical structures both at 
high-z and in local universe. 



4.4. Completeness o/H I Line Sample 

For a statistical analysis, especially the number density 
analysis, the completeness of the H I line detection is in- 
fluenced by the detection limit of absorption lines (e.g., 
equivalent width or column density). In this study, we 
have used the H I lines detected in the 40 HIRES spectra 
that have various S/N ratios. The strong line sample will 
have subtle biases arising from the selection of the quasars 
because they were once thought to be good targets for the 
detection of deuterium. For example, we avoided quasars 
with no LLSs, and we avoided LLSs with previously known 
complex velocity structure. Nevertheless, we confirmed 
that our sample is almost complete for weak lines in the 
following way. 

The minimum detectable equivalent width in the 
observed- frame, W m i n , can be estimated using the follow- 
ing relation, 

JJ = W mm N c = WrninjS/N) 

<j{W min N c ) (MlMc 1 + M L - W mm yl 2 ' 

where Ml and Mp are the numbers of pixels over which 
the equivalent width and the continuum level (Nc) are 
determined (Young et al. 1979; Tytler et al. 1987). The 
value of (S/N) is the S/N ratio per pixel. When we set 
U ~ W/a(W) is 4 (i.e., 4cr detection), the eqn.© can be 
solved to give, 

W min = (S/A)- 2 {[64+ 16(S/N) 2 x 

(Ml+mI/Mc)} 1 ^ 2 ~8}x AX (A), (4) 

where AA is the wavelength width per pixel in angstroms 
(Misawa et al. 2002). Here, we set Ml for 2.5 times 
FWHM of each line, and Mc for full width of each echcllc 
order. Once the minimum rest-frame equivalent width, 
W rest [= W m in/ (1 + z)], has been evaluated, it can be con- 
verted to the minimum column density by choosing a spe- 
cific Doppler parameter; the result is insensitive to the 
choice on the linear part of the curve of growth. Among 
the 86 H I systems in our data sample, the H I system 
at z a bs = 2.940 in the spectrum of Q0249-2212 is located 
in the region with the lowest S/N ratio (i.e., S/N ~ 11). 
This corresponds to a 4cr detection limit of log Nhi ~ 12.3 
for an isolated Lya line with any Doppler parameter seen 
in our sample (b = 15 - 80 km s _1 ). Thus, our sample 
is complete for H I lines with log Nhi > 12.3. Therefore, 
the bend in the column density distribution near log Nhi 
~ 13 is probably due to the line blending and blanketing. 

5. SUMMARY 

We present 40 high-resolutional (FWHM = 8.0 km s _1 ) 
spectra obtained with Keck+HIRES. Over the wide col- 
umn density range (12 < log Nhi < 19), we fit H I lines 
by Voigt profiles using not only Lya line but also higher 
Lyman series lines such as Ly/3 and Ly7 up to Lyman limit 
when possible. To investigate the detailed line properties, 
we made several sub-samples that are separated according 
to the distance from the quasar, redshift, the column den- 
sity, and the S/N ratio of the spectrum. We also classify 
them into HDLs (lines arising in or near to intervening 
galaxies) and LDLs (lines not obviously near to galaxies 
and hence more likely to be from the intergalactic diffuse 
gas), based on the clustering properties. The main results 
are summarized below: 
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1 . We present a database of H I absorption lines with 
a wide column density range (i.e., logNjji = 12-19) 
from a wide redshift range (i.e., z = 2-4). 

2. Our data sample is complete at log Njji > 12.3 with 
4cr line detection. The turnover at log N^i < 13 
seen in the log Njji distribution is not due to a qual- 
ity of our spectra but due to the line blending and 
blanketing. 

3. The power-law indices of the column density distri- 
bution of LDLs shows evolution with redshift, from 
P = 1.52±0.09 at z > 2.9 to f3 = 1.71±0.06 at z < 
2.9. This trend could be related to the hierarchical 
clustering in cosmological timescale. No evolution 
is seen for HDLs. 

4. Within 5,000 km s _1 of the quasars, the power-law 
index of the column density distribution for LDLs 
{fi = 1.90±0.16) is larger than those far from the 
quasars [fi = 1.53±0.05). We also found a hint (Fig- 
ure 6) that the Doppler parameters are larger near 



the quasars. These results could be due to the UV 
flux excess from the quasars. We do not see any 
similar trend for the HDLs. 

5. We suggest that HDLs and LDLs are produced by 
physically different phases or absorbers, because 
they have four key differences seen in (i) cluster- 
ing property, (ii) redshift evolution, (iii) Proximity 
effect, and (iv) \ogN HI - b min relation (see Misawa 
et al. 2004). 
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HI Lines in 40 HIRES Quasars 



Table 1 

Keck HIRES Spectra of 40 quasars 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


quasar a 


^ em 




m R 


\ min d 


\ max e 














(A) 




Q0004+1711 


2.890 


18.70 




3510 


5030 


11.9 


Q0014+8118 


3.387 




16.1 


3650 


6080 


48.8 


Q0054-2824 


3.616 




17.8 


4090 


6510 


18.2 


Q0119+1432 


2.870 


17.4 




3200 


4720 


23.7 


HE0130-4021 


3.030 


17.02 




3630 


6070 


52.5 


Q0241-0146 


4.040 


18.20 




4490 


6900 


7.5 


Q0249-2212 


3.197 


17.70 




3500 


5020 


11.0 


HE0322-3213 


3.302 


17.80 




3830 


5350 


12.8 


Q0336-0143 


3.197 




18.8 


3940 


6390 


12.7 


Q0450-1310 


2.300 


16.50 




3390 


4910 


17.1 


Q0636+6801 


3.178 




16.9 


3560 


6520 


53.4 


Q0642+4454 


3.408 




18.4 


3930 


6380 


19.0 


HS0757+5218 


3.240 


17.3 




3590 


5120 


21.5 


Q0805+0441 


2.880 


18.16 




3800 


6190 


15.7 


Q0831 + 1248 


2.734 


18.10 




3790 


6190 


29.2 


HE0940-1050 


3.080 


16.90 




3610 


6030 


35.7 




2 644 


i fi A n 




ouyu 




/IS fi 


Q1017+1055 


3.156 




17.2 


3890 


6300 


44.8 


Q1055+4611 


4.118 


17.70 




4450 


6900 


40.3 


HS1103+6416 


2.191 


15.42 




3180 


5790 


78.1 


Q1107+4847 


3.000 


16.60 




3730 


6170 


51.8 


Q1157+3143 


2.992 


17.00 




3790 


6190 


35.9 


Q1208+1011 3 


3.803 




17.2 


3730 


6170 


21.8 


Q1244+1129 


2.960 


17.70 




3370 


4880 


9.9 


Q1251+3644 


2.988 


19.00 




3790 


6190 


32.5 


Q1330+0108 


3.510 




18.56 


4030 


6450 


8.8 


Q1334-0033 


2.801 


17.30 




3730 


6170 


24.7 


Q1337+2832 


2.537 


19.30 




3170 


4710 


31.1 


Q1422+2309 h 


3.611 




15.3 


3740 


6180 


136 


Q1425+6039 


3.165 




16.0 


3730 


6170 


43.5 


Q1442+2931 


2.670 


16.20 




3740 


6180 


29.4 


Q1526+6701 


3.020 


17.20 




3460 


4980 


9.7 


Q1548+0917 


2.749 


18.00 




3730 


6180 


21.9 


Q1554+3749 


2.664 


18.19 




3240 


4770 


13.2 


HS1700+6416 


2.722 


16.13 




3730 


6180 


66.2 


Q1759+7539 


3.050 


16.50 




3580 


5050 


30.9 


Q1937-1009 


3.806 




16.7 


3890 


7450 


76.9 


HS1946+7658 


3.051 


16.20 




3890 


6300 


136 


Q2223+2024 


3.560 




18.5 


4120 


6520 


12.9 


Q2344+1228 


2.763 


17.50 




3410 


4940 


8.1 



a Quasar names are based on B1950 coordinates. 
h V magnitude from Veron-Cetty & Veron (2003). 

C R magnitude from USNO-A2.0 Catalog (Monet ct al. 1998), except for 
Q1330+0108, whose R magnitude comes from the USNO-B Catalog (Monet 
ct al. 2003). 

d Lower limit of the observed quasar spectrum. 

c Uppcr limit of the observed quasar spectrum. 

f S/N ratio at the center of each spectrum. 

s This lensed quasar is amplified by a factor of ~3.1 (Barvainis & Ivison 
2002). 

h This lensed quasar is amplified by a factor of 15.38 (Kormann ct al. 
1994). 



Table 2 



Properties of 86 H i systems 



(1) 


(2) (3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) (11) 


(12) 


(13) 


quasar 


z abs logATi a 


logiV 2 b 


NJN 2 


S/N(Ly-l) 


S/N(Ly-2) 


S/N(Ly-5) 


S/N(Ly-10) 


c d 

n iooo n S ys 


status e 


reference^ 




(cm- 2 ) 


(cm- 2 ) 



















Q0004+1711 


2 


8284 


15 


51 


14 


46 





0896 


18 


9.6 


2.6 


1.8 


10 


2 


V5000 






2 


8540 


15 


75 


14 


94 





1546 


27 


8.5 


3.7 


2.2 


17 




A. 


V5000 






2 


8707 


19 


93 


16 


03 





0001 


27 


11 


3.4 


2.0 


9 




A, 


C, V5000 


1 


Q0014+8118 


2 


7989 


18 


30 


18 


02 





5282 


63 


5.2 






11 


3 






1 




2 


9090 


16 


09 


15 


60 





3221 


88 


11 


2.0 




18 


4 






1 




3 


2277 


15 


33 


15 


28 





8855 


93 


41 


9.3 


4.2 


18 


7 






1 




q 


3212 


16 


60 


16 


24 


n 


4438 


154 


48 


14 


9 


16 


g 


V5000 




Q0054-2824 


3 


2370 


15 


56 


15 


18 





4108 


17 


8.0 






14 


7 










3 


3123 


16 


64 


14 


91 





0184 


25 


4.9 






16 


4 










3 


4488 


16 


67 


15 


21 





0346 


32 


12 


3.1 


2.2 


17 


9 










3 


5113 


15 


89 


14 


84 





0899 


50 


14 


4.9 


3.2 


17 


5 






1 




3 


5805 


17 


44 


15 


94 





0318 


96 


16 


6.3 


5.9 


21 


6 


V5000 


1 


Q0119+1432 


2 


4299 


15 


93 


15 


14 





1646 


39 


11 






6 


4 










2 


5688 


16 


39 


14 


62 





0169 


41 


16 


4.4 




11 


5 










2 


6632 


19 


37 


15 


82 





0003 


49 


25 


7.6 


5.7 


11 




C 






HE0130-4021 


2 


8581 


15 


15 


14 


84 





4902 


59 


7.2 






19 


4 








Q0241-0146 




































Q0249-2212 


2 


6745 


19 


01 


14 


28 





0000 


15 


4.8 






9 




C 




1 




2 


9401 


17 


23 


14 


65 





0026 


11 


7.8 


2.4 


3.4 


17 


4 






1 


HE0322-3213 


3 


0812 


15 


68 


14 


86 





1515 


26 


15 


6.2 




17 


3 










3 


1739 


19 


43 


14 


18 





0000 


27 


18 


11 


7.5 


9 




A, 


c 






3 


1960 


16 


61 


15 


74 





1345 


36 


22 


8.7 


6.5 


15 




A 








3 


3169 


16 


16 


15 


33 





1475 


103 


25 


13 


12 


11 




Vi 


300 





Q0336-0143 
Q0450-1310 



Q0636+6801 


2 


6825 


15 


57 


15 


02 





2837 


64 


19 






15 


3 




1 




2 


8685 


15 


85 


14 


49 





0431 


87 


41 


7.3 




10 


3 




1 




2 


9039 


18 


22 


15 


45 





0017 


60 


42 


13 


7.8 


19 


6 




1 




3 


0135 


15 


79 


14 


95 





1465 


105 


27 


17 


16 


19 


3 


D 


1 




3 


0675 


15 


28 


14 


30 





1054 


117 


44 


23 


13 


18 


6 




1 


Q0642+4454 


2 


9726 


17 


36 


14 


68 





0021 


22 


7.2 






18 


3 


D 


1 




3 


1230 


19 


48 


17 


54 





0116 


23 


8.9 






11 




C 


1 




3 


1922 


15 


27 


14 


48 





1640 


28 


15 


1.3 




18 


3 




1 




3 


2290 


15 


52 


15 


37 





7158 


27 


15 


3.4 




19 




A 






3 


2476 


16 


55 


15 


37 





0669 


29 


16 


4.0 




18 




A, B 


1 




3 


3427 


15 


40 


14 


84 





2744 


20 


17 


7.3 


4.3 


15 




B, V5000 


1 


HS0757+5218 


2 


7261 


15 


46 


14 


98 





3360 


35 


11 






10 


2 








2 


8922 


18 


34 


14 


90 





0004 


25 


18 


1.4 




13 


1 








3 


0398 


19 


82 


16 


74 





0008 


30 


25 


11 


6.7 


10 




C 




Q0805+0441 


2 


7719 


15 


30 


15 


14 





6906 


29 


5.8 






25 


7 




1 




2 


8113 


15 


99 


14 


88 





0765 


37 


8.3 






17 


4 




1 


Q0831+1248 


2 


7300 


15 


74 


14 


07 





0212 


57 


13 






11 




V1000 




HE0940-1050 


2 


8283 


10 


41 


16 


05 





4305 


52 


12 






20 


15 








2 


8610 


17 


06 


14 


53 





0029 


56 


18 


2.6 




18 


4 








2 


9174 


15 


92 


15 


35 





2669 


63 


21 


6.2 




19 


5 








3 


0387 


15 


55 


13 


84 





0196 


91 


27 


7.6 


6.7 


10 


3 


V5000 




Q1009+2956 


2 


1432 


17 


82 


15 


33 





0032 


71 


17 






9 


6 








2 


4069 


18 


80 


14 


25 





0000 


126 


36 


9.6 




9 




A 


1 




2 


4292 


17 


34 


14 


53 





0015 


109 


45 


14 


1.7 


18 




A 






2 


5037 


17 


26 


15 


49 





0167 


131 


53 


21 


15 


14 


4 




1 


Q1017+1055 


2 


9403 


15 


56 


14 


49 





0844 


12 


8.2 






11 


2 




1 




3 


0096 


15 


98 


14 


80 





0658 


35 


6.9 






21 


5 




1 




3 


0548 


17 


06 


15 


54 





0302 


25 


8.9 






18 


11 




1 




3 


1120 


15 


26 


15 


01 





5536 


43 


11 






26 


7 


Vsooo 


1 


Q1055+4611 


3 


8252 


15 


98 


15 


64 





4603 


53 


37 


15 




26 




A 






3 


8495 


16 


74 


16 


04 





1997 


31 


35 


14 


4.4 


23 




A 






3 


9343 


17 


30 


16 


32 





1035 


40 


34 


22 


13 


25 




B 





HS1103+6416 



Table 2 — Continued 



(1) 


(2) (3) (4) (5) 


(6) 


(7) 


(8) 


(9) (10) 


(11)^ (12) ^ 


(13) 


quasar 


tabs logATi a logN 2 b Ni /N 2 
(cm- 2 ) (cm- 2 ) 


S/N(Ly-l) 


S/N(Ly-2) 


S/N(Ly-5) 


S/N(Ly-10) mono c 


f^sys status 


reference^ 



Q1107+4847 


2 


7243 


16 


58 


13 


92 





0022 


38 


7.5 






12 


3 


D 








2 


7629 


19 


13 


17 


51 





0239 


43 


12 






12 




C 


1 






2 


8703 


15 


25 


14 


76 





3226 


83 


18 






17 


7 








Q1157+3143 


2 


7710 


17 


63 


14 


50 





0009 


69 


22 






13 


3 










2 


8757 


15 


66 


15 


54 





7713 


85 


28 






19 


9 










2 


9437 


17 


44 


17 


16 





5282 


99 


40 






18 


5 


V5000 






Q1208+1011 


3 


3846 


17 


35 


15 


04 





0049 


24 


18 


6.1 


3.0 


19 


6 










3 


4596 


16 


88 


16 


03 





1430 


22 


18 


8.4 


5.5 


19 


10 










3 


5195 


10 


15 


15 


73 





3802 


26 


19 


11 


7.7 


24 


6 










3 


7206 


15 


48 


14 


65 





1485 


27 


19 


14 


12 


21 


6 








Q1244+1129 


2 


9318 


15 


97 


14 


87 





0784 


31 


14 


5.0 


3.5 


17 


3 


V5000 






Q1251+3644 


2 


8684 


15 


82 


14 


03 





0161 


63 


20 






14 


3 








Q1330+0108 




































Q1334-0033 


2 


7572 


15 


40 


14 


24 





0693 


60 


9.0 






15 


3 


V5000 






Q1337+2832 


2 


4336 


18 


92 


16 


32 





0025 


60 


14 


2.5 




14 


8 










2 


5228 


15 


81 


14 


44 





0423 


161 


22 


6.9 


2.1 


18 


5 


V5000 






Q1422+2309 


3 


3825 


16 


53 


16 


33 





6427 


389 


278 


83 


54 


19 


4 










3 


5362 


15 


94 


15 


83 





7691 


462 


370 


170 


70 


22 




B, V5000 






Q1425+6039 


2 


7700 


19 


37 


16 


20 





0007 


38 


6.3 






17 




C, D 




1 1 




2 


8258 


20 


00 


19 


68 





4716 


33 


9.0 






7 




C 


1 






3 


0671 


16 


20 


14 


90 





0496 


94 


13 


4.4 


1.1 


16 


3 










3 


1356 


10 


66 


16 


15 





3107 


184 


127 


6.4 


2.7 


17 




B, V5000 




CD 


Q1442+2931 


































CZ3 


Q1526+6701 


2 


9751 


15 


12 


15 


11 





9808 


24 


7.5 


4.5 


3.0 


13 


5 


V5000 




B' 


Q1548+0917 




































Q1554+3749 


2 


6127 


17 


97 


14 


45 





0003 


43 


7.5 


2.7 


1.5 


11 


3 


V5000 







HS1700+6416 


































K 


Q1759+7529 


2 


7953 


15 


26 


14 


92 





4584 


44 


22 






16 


5 




1 


1 — 




2 


8493 


17 


44 


15 


60 





0145 


49 


27 


7.9 




16 


6 




1 






2 


9105 


19 


90 


17 


62 





0052 


60 


32 


10 


6.8 


15 




C 


1 


LTj 


Q1937-1009 


3 


5725 


17 


94 


15 


89 





0089 


147 


61 


27 


20 


20 


4 




1 


HS1946+7658 


3 


0498 


17 


45 


14 


42 





0009 


266 


44 






10 




D. Vionn 


1 


.O 


Q2223+2024 




































Q2344+2024 


2 


4261 


18 


46 


15 


18 





0005 


16 


6.8 






8 


3 




1 


9= 




2 


6356 


15 


45 


14 


03 





0382 


21 


11 


1.5 




7 


2 




1 


-A 




2 


7107 


10 


64 


15 


68 





1092 


31 


13 


4.8 


1.7 


12 


6 


V5000 


1 






2 


7469 


16 


67 


16 


27 





4039 


34 


14 


6.7 


3.1 


19 


10 


V5000 







a H I The largest column density. 

b H I The second largest column density within ±1000 km s _1 of the largest column density. 
c Numbcr of H I lines within ±1000 km s _1 of the main component. 
d Numbcr of HDLs (sec § 3.3). 

C A : centers of H I systems are separated by less that 1000 km s _1 , B : there are gaps in the spectrum within 1000 km s~ 1 of the main 
component, C : normalization of spectrum is not good because of strong absorption with log Nhi > 19, D : candidate for a quasar intrinsic 
system (Misawa et al. 2007), Viooo : velocity difference from the quasar emission redshift is smaller than 1000 km s _1 , V5000 : velocity 
difference from the quasar emission redshift is smaller than 5000 km s _1 . 

f 1 : Listed in NED or literature in 2002 (Peroux et al. 2001; Storrie-Lombardi et al. 1996; Petitjcan, Rauch, & Carswell 1994; Lu ct al. 
1993; Steidel & Sargent 1992; Lanzetta et al. 1991; Barthcl, Tytler, & Thomson 1990; Stcidcl 1990a,b; SBS; SSB). 



Misawa et al. 



Table 3 

Sub-samples of H I lines for statistical analysis 



(1) 


(2) 


(3) 


(4) 


Sub-samplc a 


Nays b 


N HnB c 


Criteria 


SO 


86 


1339 


All H I systems and lines 


SI 


61 


973 


H I systems meeting the conditions in § 2 d 


S2a 


48 


767 


Av(z cm - z ab3 ) > 5000 km s -1 


S2b 


13 


206 


Au(z cm - z abs ) < 5000 km s _1 


S3a 


34 


554 


S/N ratio at Lya is larger than 40 


S3b 


17 


280 


S/N ratio at Lya is larger than 70 


S4a 


30 


419 


Zab S < 2.9 


S4b 


31 


554 


Zabs > 2.9 


^51213 




244 


12 < log Nhi < 13 


S5l214 




716 


12 < log N HI < 14 


S5i215 




806 


12 < logN HI < 15 


S5l216 




885 


12 < logN HI < 16 


S5i3ig 




728 


13 < logN HI < 19 


S5i419 




256 


14 < log Nhi < 19 


S5i5ig 




106 


15 < log N H i < 19 


S5iei9 




58 


16 < log Nhi < 19 



a Sub-samplcs S2a — S5iei9 are all derived from a sample SI. 

b Number of H I systems, i.e. lines with log Nhi > 15 cm -2 . 

c Number of H I lines, i.e. all H I lines within ±1000 km s~ x of the systems. 

d Should satisfy four conditions; (i) separated from the quasar by more than 
1000 km s _1 , (ii) column density of the main component is logNni < 19 cm -2 , (iii) 
there arc no gaps in the HIRES spectrum near Lya, and (iv) there are no other H I 
systems within ±2000 km s _1 . 



Table 4 



Parameters of the column density distribution 



(1) 

sub-sample 


(2) 
a 


(3) 
A b 


(4) 
Az c 


SI (all lines) 


1.398 ± 0.025 


7.389 ± 0.407 


1.606 


S2a (Ati > 5000 km/s) 
S2b (Av < 5000 km/s) 


1.390 ± 0.027 
1.360 ± 0.048 


7.262 ± 0.431 
6.801 ± 0.753 


1.266 
0.341 


S3a (S/N > 40) 
S3b (S/N > 70) 


1.368 ± 0.034 
1.350 ± 0.040 


6.937 ± 0.552 
6.660 ± 0.615 


0.887 
0.453 


S4a (z < 2.9) 
S4b (z > 2.9) 


1.343 ± 0.031 
1.439 ± 0.031 


6.537 ± 0.494 
8.015 ± 0.497 


0.741 
0.865 


Kirkman & Tytlcr (1997) 
Pctitjean et al. (1993) 


1.5 

1.46 


8.79 
8.08 





a Bcst fit value and la error of (3 in cqn.Q. 
b Bcst fit value and ltr error of A in can. fl} . 
c Total rcdshift width of sub-sample. 



Table 5 
k-s test for column density 



(1) 


(2) 


(3) 


sub-samples 


D 


Prob" 




(%) 


S2a / S2b 


0.063 


52.6 


S3a / S3b 


0.073 


21.1 


S4a / S4b 


0.055 


44.5 



a Probability that the two distri- 
butions were drawn from the same 
parent population. 
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Table 6 

k-s test for doppler parameter 



(1) 


(2) 


(3) 


sub-samples 


D 


Prob a 




(%) 


S2a / S2b 


0.116 


2.4 


S3a / S3b 


0.036 


97.1 


S4a / S4b 


0.063 


28.9 



a Probability that the two distri- 
butions were drawn from the same 
parent population. 



Table 7 



Sub-samples separated into HDLs and LDLs 





(1) 


(2) 


(3) 




(4) 






Sub-sample 


N sys a 


N Hne b 




Criteria 






SI 


HDLs 


61 


306 


H I systems 


meeting the conditions 


in 


§ 4.1 




LDLs 


61 


667 


H I systems meeting the conditions 


in 


§ 4.1 


S2a 


HDLs 


48 


240 


Av(z e . 


in- 2ab,) > 5000 km s~ 


-1 






LDLs 


48 


527 


Av(z e 


in- tabs) > 5000 km s" 


-1 




S2b 


HDLs 


13 


66 


Ati(z e 


m— Zabs) < 5000 km s~ 


-l 






LDLs 


13 


140 


Ad(z c 


iB- Zabs) < 5000 km s~ 


-l 




S4a 


HDLs 


30 


143 




z abs < 2.9 








LDLs 


30 


276 




z abs < 2.9 






S4b 


HDLs 


31 


163 




z abB > 2.9 








LDLs 


31 


391 




Zabs > 2.9 







a Number of H I systems 
b Number of H I lines 



Table 8 

Parameters of column density distribution 



(1) 




(2) 


(3) 


(4) 


sub-sample 




a 


A b 


Az c 


SI 


HDLs 


1.269 ± 0.034 


5.736 ± 0.545 


0.415 


(all lines) 


LDLs 


1.589 ± 0.075 


10.04 ± 1.045 


1.192 


S2a 


HDLs 


1.264 ± 0.032 


5.663 ± 0.512 


0.332 


(An > 5000 km/s) 


LDLs 


1.526 ± 0.054 


9.172 ± 0.757 


0.933 


S2b 


HDLs 


1.167 ± 0.047 


4.203 ± 0.739 


0.082 


(Ac < 5000 km/s) 


LDLs 


1.897 ± 0.163 


14.25 ± 2.278 


0.259 


S4a 


HDLs 


1.223 ± 0.031 


5.050 ± 0.496 


0.190 


(z < 2.9) 


LDLs 


1.712 ± 0.055 


11.64 ± 0.766 


0.551 


S4b 


HDLs 


1.283 ± 0.044 


5.933 ± 0.693 


0.227 


(z > 2.9) 


LDLs 


1.517 ± 0.091 


9.102 ± 1.274 


0.638 


Pctitjean ct al. (1993) 


(HDLs + LDLs) 


1.46 


8.08 




Kirkman & Tytler (1997) 


(LDLs) 


1.5 


8.79 





a Bcst fit value and la error of (3 in eqn.JTJ. 
b Bcst fit value and la error of A in eqn. (|T)l ■ 
c Total rcdshift width of sub-sample. 
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Table 9 

k-s test for doppler parameter 





(1) 

sub-samples 




(2) 
D 


(3) 
Prob a 

(%) 




HDLs (SI) / LDLs (SI) 




0.049 


67.9 


HDLs (Av > 5000 kms" 1 ) / LDLs (Av > 5000 km s" 
HDLs (Av < 5000 km s _1 ) / LDLs (Av < 5000 km s" 


- 1 ) 


0.084 
0.158 


18.8 
19.4 


HDLs (Av 
LDLs (Av 


> 5000 km s" 1 ) / HDLs (Av < 5000 km s~ 

> 5000 km s" 1 ) / LDLs (Av < 5000 km s" 


- 1 ) 
- 1 ) 


0.091 
0.143 


76.6 

2.0 




HDLs (z < 2.9) / LDLs (z < 2.9) 
HDLs (z > 2.9) / LDLs (z > 2.9) 




0.095 
0.086 


34.3 
34.6 




HDLs (z < 2.9) / HDLs (z > 2.9) 
LDLs (z < 2.9) / LDLs (z > 2.9) 




0.134 
0.041 


11.7 
94.7 



Probability that the two distributions were drawn from the same parent pop- 
ulation. 
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Fig. 1. — The solid and dashed histograms represent the numbers of H I systems in samples SO and SI as a function of redshift. The 
dot-dashed histogram is the number of quasars in which absorption lines could have been detected (namely the number of spectrum windows). 
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12 14 16 18 



log Nttt (cm 2 ) 

Fig. 2. — Observed column density distribution per unit redshift and unit column density for sample SI. The H I lines have been binned 
into intervals of 0.5 in log TV fjj. The open squares and vertical bars represent the observed data and la errors. The solid line is the best fit 
power law for our study. Dashed and dotted lines are the best fit power laws in the range of 12 < log Nhi < 14 (KT97) and 12 < log Nhi 
< 22 (Petitjean ct al. 1993). 




20 40 60 80 100 ^ 20 40 60 80 100 



b (km s 1 ) b (km s 1 ) 

Fig. 3. — Observed Dopplcr parameter distributions for sub-samples SI and S1213. The H I lines have been binned into intervals of Ab = 
5 km s _1 . The open squares and vertical bars represent the observed data with their la errors. Dashed and dotted lines are the input data 
with truncated Gaussian distribution profiles for H95 (60 = 28 km s _1 , at = 10 km s _1 , and b C ut = 20 km s _1 ) and L96 (bo = 23 km s _1 , 
(T;, = 8 km s _1 , and b cu t = 15 km s _1 ). 
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Fig. 4. — Two point correlation functions for sub-samples S5i2i5 and £5i5ig. The bin size is 50 km s _1 . Solid histogram and vertical bars 
in each bin represent the value of correlation degree, £ (v), and the Poisson error. Dotted vertical line denotes the velocity separation at which 
the lower l<r deviation of £(v) first goes below £(v) = over v > 50 km s _1 . 




00 



b (km s ) b (km s ) 

Fig. 5. — Observed Doppler parameter distributions of LDLs for sub-samples S2a (Au > 5000 km s -1 ) and S2b (Ad < 5000 km s -1 ). The 
H I lines have been binned into intervals of Ab = 5 km s -1 . Open square and vertical bars represent observed data and the la errors. Dashed 
and dotted lines are the input data with truncated Gaussian distribution profiles for H95 (bo = 28 km s _1 , cr;, = 10 km s _1 , and b m i„ = 
20 km s _1 ) and L96 (6q = 23 km s _1 , a^, = 8 km s _1 , and b m i„ = 15 km s _1 ). 
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b (km/s) 

Fig. 6. — The cumulative distribution functions for the b-values of LDLs. The 527 LDLs from Av > 5000 km s~ 1 (sample S2a) are shown 
with the dashed line, the 140 LDLs at Av < 5000 km s~ x (S2b) with the dotted line, and we show all 667 LDLs (SI) with the central thin 
solid line. Ly« forest lines near to a quasar tend to have larger Dopplcr parameters compared with lines far from a quasar. 
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log N m (cm 2 ) 

Fig. 7. — Column density distributions of LDLs at z < 2.9 (S4a) with open squares and solid lines, and at z > 2.9 (S4b) with filled triangles 
and dashed lines. The line frequency of stronger LDLs (i.e., logNjfj > 14.5) at z < 2.9 preferentially decreases compared with those of LDLs 
at z > 2.9, while the frequency of weaker LDLs does not change with redshift. 
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HDLs 
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Fig. 8. — Cartoon of the distribution of absorbers at z > 2.9 and z < 2.9. Shaded circles are HDL absorbers. Dotted and small circles are 
strong (e.g., log Nhi > 14.5) and weak (e.g., log Nhi < 14.5) LDL absorbers, respectively. If strong LDL absorbers would gather around HDL 
absorbers within relative velocity of Av < 200 km s _1 , the number of strong LDLs decrease as redshift decreases. Such trend is consistent 
with the concept of the hierarchical clustering model. 
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log N ffl (cm 2 ) 

Fig. 9. — Column density distributions of LDLs far from the quasars with relative radial velocity of Ad > 5000 km s _1 (S2a) which are 
shown as open squares and solid lines, and near the quasars with Ad < 5000 km s _1 (S2b) denoted as filled triangles and dashed lines. 
The line frequency of stronger LDLs near the quasars preferentially decreases compared with those of LDLs far from the quasars, while the 
frequency of weaker LDLs is not affected by the distance from the quasars. 



hot region 




Fig. 10. — Cartoon of the distribution of absorbers far from the quasars (Ad > 5000 km s" 1 ) and near the quasars (Ad < 5000 km s" 1 ). 
Shaded regions are dense regions that arc adiabatically compressed and heated, while dotted regions arc diffuse cold regions. If LDL absorbers 
near the quasars are strongly affected by the UV flux from the quasars, H I gas at the diffuse cold regions would be preferentially ionized. As 
a result, stronger LDLs would become to be weaker LDLs and only central hot regions are observed, which makes the relative frequency of 
stronger LDLs smaller and the mean b value of LDLs larger. 
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APPENDIX 
DISCUSSION OF INDIVIDUAL H I SYSTEMS 

In this section, we describe the results of fitting the 86 H I systems in sample SO. Velocity plots of them with ± 
1000 km s _1 widths for the lowest five orders of Lyman series (i.e., Lya, Ly/3, Lyy, LyS, and Lye) are presented in 
Figure 11 as far as they are accessible. In Table 10 we give in column (1) ID number; columns (2) and (3) observed 
wavelength and velocity shift from the system center; column (4) absorption redshift; columns (5) and (6) column density 
with la error; columns (7) and (8) Doppler parameter with la error; column (9) line identification. If narrow lines with 
b < 15 km s _1 are not identified as specific metal lines, we use "M i" as unidentified lines in the column (10). Table 10 
lists only H i, M I, and metal lines that are detected within ±1000 km s _1 windows of the 81 H I systems. Metal lines in 
the H I system windows are neither numerated in the table nor marked with ticks in Figure 11 because they happen to 
locate within the H I system windows and they are not physically relate to the H I systems. Important metal absorption 
lines in the 86 H I systems that are detected in our spectra are also summarized in a separate table (Table 11). 

Q0004+1711 (z em = 2.890). - SSB observed this quasar, and detected strong C iv and Si iv absorption lines at 
z a 6;5=2.5181 as well as a strong Mg n line at z a b s = 0.8068. We confirm the prominent LLS at z a bs — 2.881. We 
see Si n A1260, Si n A1527, C n A1335, and O I A1302 lines, but no C iv doublet. Our spectrum has a range of 
3510 A to 5030 A. Both Lya and Ly/3 are detected at z abs = 2.422 - 2.890. 

z a bs = 2.8284 — Although the spectrum has a range of Lya up to Lyl3, the S/N ratio is very low (S/N = 18 at 
Lya, and 1.8 at LylO). This system is within 5000 km s _1 of the emission redshift of the quasar at z em = 2.89. 

Zabs = 2.854-0 — This system is also within 5000 km s _1 of the quasar. Though the spectrum covers the Lyman 
limit of the system (\u m it ~ 3513 A), the low S/N ratio of the spectrum prevented us from measuring this. This 
system is shifted only 1300 km s^ 1 blueward of the DLA system at z a b s = 2.8707, and Lya is strongly blended 
with the left wing of the DLA. 

Zabs = 2.8707 — This system was previously detected by SBS. Most components in the system are blanketed by 
the wings of the main component, which has a large column density, log Nhi — 19.93, and rather small Doppler 
parameter, b = 12.57 km s _1 . This system is also within 5000 km s _1 of the quasar. 

Q0014+8118 (z em = 3.387). — This quasar has been well studied since its discovery in 1983 (Kuhr et al. 1983), as 
there is a candidate D 1 line at z a t, s — 3.32. An upper limit on the D/H ratio was determined to be D/H < 25 - 
60 xl0~ 5 for this system (Songaila et al. 1994; Carswell ct al. 1994). Rugen & Hogan (1996a, b) also detected a 
D 1 line in another LLS at z a bs = 2.80 in this quasar. Buries, Kirkman, & Tytler (1999), however, claimed that 
these absorption lines were not primarily due to D 1, based on their improved spectrum. Our spectrum ranges from 
3650 A to 6080 A. Both Lya and Ly/3 are detected at z abs = 2.558 - 3.387. 

Zabs = 2.7989 — The absorption profile around the main component (the "central trough" hereafter) is strongly 
damped for Lya, Ly/3, and Ly7, which makes it difficult to fit the profile. If the trough was fit with a single 
component, the Doppler parameter was found to be rather large, b > 60 km s _1 . Fortunately this system has many 
C iv and Si iv lines. Therefore the C iv lines were used as a reference, and the trough was fit with two components 
having b = 45 and 33 km s _1 , respectively. Both components are found to have high column densities, log Nhi > 
18. They may be resolved into narrower components. 

Zabs = 6 2- 9090 — If the central trough was fit with only one line, the column density was found to be log Nhi > 
16.8. However, there is no Lyman break feature around 3565 A. Ly/3 has an asymmetrical profile. Therefore we 
fit the trough with two components having log Nhi = 16.09 and 15.60. The best-fitting model for Lya and Ly/3 is 
slightly inconsistent with Lyy and LyS. 

Zabs = 3.2277 — This is a very weak system with log Nhi = 15.33, that may be a strong Lya forest member 
produced by an intergalactic cloud. There are no metal lines in the system. The spectrum has a narrow data defect 
at Av = —500 to —400 km s~ 1 from the main component in Lya window. 

Zabs — 3.3212 — Buries, Kirkman, & Tytler (1999) fit the central trough with four components positioned at 
Aw = —98.6, 0, +98, and +155 km s~ x from the main component. We also fit the trough with four components 
positioned at Av =—90, 0, +100, and +150 km s _1 from the main component, which is in good agreement with 
the results of Buries, Kirkman, & Tytler (1999). The Lyman break around 3940 A suggests that this system has a 
column density larger than log Nhi > 16.6. This system is within 5000 km s _1 of the quasar. A C iv complex at 
z a bs = 2.40 around 5260 - 5275 A is blended with Lya lines of this system. 

Q0054-2824 (z em = 3.616). — In this quasar, SSB found strong Mg 11 lines at z abs = 1.3412 and 1.4398, and three Si iv 
lines at z a bs = 3.2791, 3.5068 and 3.5800, but associated C iv lines were not detected. The Si iv system at z a bs = 
3.5800 is known to be associated with the conspicuous LLS at z a bs = 3.585. Our spectrum ranges from 4090 A to 
6510 A. Both Lya and Ly/3 are detected at z a ts — 2.987 - 3.616. 

Zabs = 3.2370 — This is a less reliable system, because the spectrum contains only three Lyman lines (Lya, Ly/3, 
and Ly7), and the S/N ratio is very low (S/N = 17 at Lya). If the central trough is fit with only one component, 
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the column density of the main component is log Njji = 16.6. The ratio of Ni (the largest H I column density 
in the fit) to N2 (the second largest H 1 column density in the fit) is then N1/N2 > 200, which is rather large 
compared with the usual value of Ni/N 2 ~ 0.3. Therefore we used two components to fit the trough. We detected 
Si in and Si iv in this system. 

z a bs = 3.3123 — The model fit for Lya, Ly/3, and Ly7 is inconsistent with Ly5, which may be due to the low 
S/N ratio of the spectrum around the Ly<5 lines. There are no corresponding metal lines, despite the large column 
density of the main component, log N hi > 16. 

z a bs = 3.4488 — H I lines of orders higher than LyS are located in the Lyman continuum of the LLS at z a t, s = 
3.58. There are no corresponding metal lines at this redshift. 

Zabs = 3.5113 — SSB detected a tentative Si iv doublet at z a bs — 3.507. Referring to that line, we found a H 1 
system with log N hi = 15.9 at z a bs — 3.511. There are seven unidentified narrow lines with 7 < b < 14 km s" 1 at 
5487 - 5499 A. 

Zabs = 3.5805 — This LLS was reported in SBS. Since many components are heavily blended with each other 
in the central trough, the profile of higher order lines were used as a reference, and the profile was fit with four 
components. The Lyman break of this system is detected around 4180 A. This system has various metal lines, 
including C 11, Si 11, and Si iv. We did not find any velocity shift between low-ionization ions (C 11 and Si 11) 
and high-ionization ions (Si iv), though such velocity shifts are often seen in DLA systems (e.g., Lu et al. 1996b; 
Prochaska et al. 2001). This system is within 5000 km s _1 of the quasar. 

Q0119+1432 (z em — 2.870). — This quasar was discovered during the course of the Hamburg/CfA Bright Quasar 
Survey (Hagen et al. 1995; Dobrzycki et al. 1996). No detailed analysis of this quasar has been published. Our 
spectrum ranges from 4090 A to 6510 A. Both Lya and Ly/3 are detected at z a bs = 2.120 - 2.870. 

z a bs = 2.4299 — Since this system is at low redshift compared to other H 1 systems in this study, the number density 
of Lya forest lines is relatively small around this system. The spectrum covers only Lya and Ly/3. Nonetheless 
the model is reliable thanks to the low number density of Lya forest lines. Only six H 1 lines are detected within 
1000 km s _1 of the main component. 

Zabs = 2. 5688 — The central trough of this system has a simple profile, but the profiles of higher-order lines suggest 
that this system has a complex structure of narrow H 1 components. We used four components to fit the trough. 
We did not detect any metal lines in the system. 

Zabs = 2.6632 — The large column density of the main component, log N hi — 19.37, results in strong Doppler 
wings at both sides of the line. A Lyman-break feature is also detected around 3350 A. We detected low-ionization 
Si 11 and Si 111 lines, while high-ionization lines such as Si iv and C iv were not detected. 

HE0130-4021 (z em = 3.030). — This quasar was discovered by Osmer & Smith (1976). Kirkman et al. (2000) obtained 
the spectrum of the quasar with a total integration time of 22 hr, and found an LLS at z aos — 2.8 with low D/H 
abundance ratio, D/H — 3.4 x 10~ 5 . The spectrum covers the range 3630 A to 6070 A. Both Lya and Ly/3 can be 
detected at z a b s = 2.539 - 3.030. 

Zabs = 2.8581 — Though the spectrum covers Lya, Ly/3, Lyy, and Ly5, they are located in a region with low 
S/N ratio, S/N < 10. The main component has a column density just over the limiting value of logNni = 15. 
Nonetheless, this system is probably not a normal Lya forest line, because it is accompanied by many metal lines 
such as C iv, Si iv, Si 111, and Si n. 

Q0241-0146 (z em = 4-040). — The emission lines of the quasar, such as Lya, O 1, C 11, Si iv, and C iv , are known 
to be very broad and rounded. Storrie-Lombardi et al. (1996) found a DLA system at z a bs — 2.86 with logNni 
= 19.8 and a metal line system with Mg 11 and Fe 11 lines at z aos = 1.435. Our spectrum ranges from 4490 A to 
6900 A. Both Lya and Ly/3 were detected at z aos = 3.377 - 4.040. However, no H 1 lines with logNni > 15 were 
detected in our spectrum. 

Q0249-2212 (z em = 3.197). — SSB found a very strong C iv system at z a bs =3.1036 and weak C iv systems at z aos 
=2.6736 and 3.1294. The LLS detected at z a bs — 2.937 has no associated metal lines. Our spectrum ranges from 
3500 A to 5020 A. Both Lya and Ly/3 are detected at z abs = 2.412 - 3.129. 

z a bs = 2.6745 — This system is a sub-DLA system with a column density logNni = 19.0, resulting in strong 
Doppler wings; however, the S/N ratio is very low (S/N = 15 at Lya). The Lyman break of this system at A = 
3349 A is not covered by our spectrum. 

Zabs = 2.9401 — This system is also detected in a region of low S/N ratio (S/N = 11 at Lya). The main component 
has logNni = 17.2. There is a tentative Lyman break feature around 3595 A. This system, however, does not have 
any metal lines, in spite of the large column density; this has already been noted by SSB. 

HE0322-3213 (z em = 3.302). — Coordinate of this quasar is given in Kirkman et al. (2005). Our spectrum ranges from 
3830 A to 5350 A. Both Lya and Ly/3 are detected at z abs = 2.734 - 3.317. 
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z a bs — 3.0812 — We fit the central trough with one component having log N hi — 15.68. The lines around the 
main component have column densities similar to that of the main component; log Nhi = 14.81, 14.19, 14.48, and 
14.86 at Av = —600, +250, +550, and +900 km s _1 from the main component. 

z a bs = 3.1739 — Our spectrum covers from Lya to Lyll of this system. It was not possible to fit the central 
trough well, as it is asymmetrical. This effect is probably artificial, because it seems to be caused by the failure of 
continuum fitting, as is often the case for the spectrum around strong absorption features. Nonetheless, our fitting 
model is reliable to some extent, since the profiles of higher orders are fit very well. This system is accompanied 
by five Si n lines. 

z a bs — 3.1960 — This system is only 1600 km s _1 redward of the system at z aos — 3.1739. Various metal lines, 
such as C ii, Si n, and Si in, were detected in the system. Though the spectrum has a wide data defect between 
5126 A and 5132 A, corresponding to Av — 600 - 900 km s _1 from the main component, we were able to fit the 
regions with reference to the profiles of higher orders. 

z a bs = 3.3169 — Our spectrum covers from Lya to Lyman limit of this system, and the S/N ratio is very high 
(e.g., S/N = 103 at Lya). There are many broad and smooth components blueward of the main component, while 
only narrow lines were detected in regions redder than the main component. The narrow line clustering around 
5262 A corresponds to unidentified metal lines. This system has corresponding C n and C in lines. 

Q0336-0143 (z em = 3.197). — This quasar was discovered in the course of the Large Bright Quasar Survey (LBQS), 
and is known to have a DLA system at z abs — 3.061 with log Nhi = 21.18 (Lu et al. 1993). Fe n, Si n, Si in, Si iv, 
C ii, Al ii, and O i are associated with this DLA system. Lu et al. (1993) also detected an Na i doublet at z a bs = 
0.1666, though this identification is less reliable because of line blending with Al n A1671 at z abs = 3.1146. The 
Mg n doublet at z a bs — 1-456 is also uncertain, because the blue and red members of the doublet poorly agree in 
rcdshift. The system also provides accurate measurements of uncommon metal lines such as Ar i, P n and Ni n 
(Prochaska et al. 2001). Our spectrum ranges from 3940 A to 6390 A. Both Lya and Ly/3 are detected at z a bs 
= 2.841 - 3.197. Unfortunately, the low S/N ratio of the spectrum prevents us from detecting not only this DLA 
system but also other H I systems with log Nhi > 15. 

Q0450-1310 (z em = 2.300). - This quasar was discovered by C. Hazard, and first studied by SBS and Steidel & 
Sargent (1992). They found one Fc 11 line at z abs = 1.1745, three Mg 11 lines at z abs = 0.4940, 1.2291 and 1.3108, 
and three C iv and Si iv lines at z a bs — 2.0669, 2.1063 and 2.2315. Petitjean et al. (1994) found an additional 
Mg 11 line at z ab s = 0.548, and four C iv lines at z abs = 1.4422, 1.5223, 1.6967 and 1.9985. The system at z abs 
= 2.0669 is a DLA candidate, because it has a strong Lya line with large rest-frame equivalent width (W res t = 6 
A), and corresponding O 1 line which is often detected in DLA systems. The system at z a bs — 2.2315 is probably 
associated with the quasar, as the velocity difference between the two is only 2080 km s , and because the system 
has a high-ionization N v doublet which is usually detected in the systems physically associated to quasars. Our 
spectrum ranges from 3390 A to 4910 A, and most of the region is redder than the peak of the Lya emission lines at 
the rcdshift of the quasar, z a bs = 2.300. Therefore we could not detect any H 1 lines in our spectrum with log Nhi 
> 15. 

Q0636+6801 (z em =3.178). — This quasar at z em = 3.178 is one of the most luminous quasars known, and is listed as 
a radio quasar in Hewitt & Burbidge (1987). There are several C iv absorption systems at z a bs — 2.4754, 2.8051, 
2.9040, 3.0174, and 3.0589. The system at z abs = 2.9040 is associated with the LLS at z ab s = 2.909. At lower 
redshift, the Mg n system is detected at z a bs = 1-2941 (SSB). Our spectrum ranges from 3560 A to 6520 A. Both 
Lya and Ly/3 were detected at z a bs = 2.471 - 3.178. 

z a bs = 2.6825 — The absorption lines around this system were well fit due to the high S/N ratio (S/N = 64 at 
Lya) and the low number density of Lya forest lines around the system. The Ly7 lines of the system are blanketed 
by the Lyman continuum of the LLS at z a bs =2.904. 

Zabs = 2.8685 — The spectrum includes the Lya to Ly8 lines of the system, though Ly7 and Ly8 are blanketed 
by the Lyman continuum of the LLS at z abs = 2.904. Four lines between Av = 400 km s _1 and 700 km s _1 from 
the main component in Lya window are probably not H I lines, since corresponding lines of higher orders are not 
detected. 

Zabs — 2. 9039 — This system is an LLS with a large column density, and has a clear Lyman break around 3570 A. 
Songaila & Cowie (1996) have already evaluated the column density of the system, finding log Nhi — 17.8. Our 
best fit to this line gave log Nhi = 18.22. We detected 8 C iv and 6 Si iv doublets in the system, though Si iv A1394 
components are affected by the spectrum gap. The system also has three O 1 lines, which strongly suggests that it 
is in a low ionized state, surrounded by gas clouds of large column density. 

Zabs = 3.0135 — This system, with \ogNni = 15.79, has two C iv doublets. The spectrum ranges from Lya to 
the Lyman limit of the system. Misawa et al. (2007) identified this system as a quasar intrinsic system, based on 
the partial coverage analysis of the C iv doublet. 
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z a bs = 3.0675 — This is a weak system with log Nhi = 15.28. The fitting model is very reliable because the 
spectrum ranges from Lya to the Lyman limit of the system with a high S/N ratio (e.g., S/N = 117 at Lya, 44 at 
Ly/3, and 13 at LylO). 

Q0642+4454 (zem—3.408). — This quasar is one of the quasars for which LLS absorption was detected for the first time 
(Carswell et al. 1975). However, we did not detect the LLS which had been discovered at z aos = 3.295 by Carswell 
et al. (1975). SSB found three C iv systems at z a bs — 2.9724, 3.1238, and 3.2483, and one Mg 11 system at z aos = 
1.2464. Our spectrum ranges from 3930 A to 6380 A. Both Lya and Ly/3 are detected at z aos = 2.831 - 3.408. 

Zabs — 2.9726 — The model is less reliable here, because only Lya and Ly/3 could be used as reference lines. 
Nonetheless, this system would be expected to have a large column density, as there are various metal lines such 
as C iv, Si 11, and Si iv. We fit the central trough with one component having log Nhi = 17.36; however, the ratio 
of iVi to N 2 is too large, Ni/N 2 > 10 3 . This component may be resolved into multiple narrow components. This 
system was classified into a quasar intrinsic system (Misawa et al. 2007). 

Zabs = 3.1230 — This system is a sub-DLA with log Nhi — 19.48. The existence of O 1 line in the system strongly 
suggests that the system has large column density, because O 1 lines are rarely seen except in DLAs. The Lyman 
break of the DLA system is seen around 3750 A in the low resolution spectrum in SSB. We did not detect C iv 
lines in this system, though Si iv and C 11 lines were detected. 

Zabs = 3.1922 — This weak system does not have any corresponding metal lines. Around Av = 700 km s _1 from 
the main component, there are two unidentified metal lines. The sharp spike at Av — —150 km s _1 from the main 
component is a data defect. 

Zabs = 3.2290 — There are two components with very similar column densities, log Nhi = 15.52 and 15.37, at 
5145 A and 5150 A. The ratio of N± to N 2 is near unity. This system is only 1300 km s _1 blueward of the system 
at Zabs = 3.248. 

Zabs = 3.2476 — Based on the C iv doublet at z aos = 3.248 detected by SSB, we found a corresponding H 1 line with 
log Nhi > 15. Though the central trough is damaged by a wide data defect of 2 A width, it was possible to fit the 
trough using the features of higher orders. The most interesting feature of this system is that the high-ionization 
lines (e.g., Si iv) are surrounded by the low-ionization lines (e.g., C 11 and C 111), which is the reverse of the trend 
usually seen in DLA systems. 

Zabs = 3.3427 — Since the Lya lines of this system are also affected by the 2 A-wide data defect, we fit the system 
by referring to the profiles of higher orders. The system is within 5000 km s _1 of the quasar. 

HS0757+5218 (z em = 3.240). — This quasar was discovered during the course of the Hamburg/CfA Bright Quasar 
Survey. No detailed analysis of the quasar spectrum has been published. Our spectrum ranges from 3590 A to 
5120 A. Both Lya and Ly/3 are detected at z abs = 2.500 - 3.212. 

Zabs = 2.7261 — Our spectrum covers the Lya, Ly/3, and Ly7 lines of the system, though Ly7 is not useful because 
of low S/N ratio. Four narrow components were detected at 600 km s _1 - 800 km s _1 from Lya of the main 
component. They are unidentified metal lines. 

Zabs = 2.8922 — The central trough has a wide and smooth profile. We fit it with a single component having 
log Nhi = 18.34. However, the Doppler parameter of the component is too large, b = 54 km s _1 . This component 
may be resolved into multiple narrow components. 

Zabs = 3.0398 — Since this sub-DLA system has strong damping wings, the spectrum could not be normalized 
correctly around the Lya line, as is often the case for echelle-formatted spectra. Nonetheless, we fit the strongly 
damped feature with five components by referring to higher orders. Metal lines corresponding to this system are 
not detected, despite the large column density of the system. 

Q0805+ 0441 (z em = 2.880). — This quasar is well known as the radio source 4C 05.34. Chen et al. (1981) first studied 
the absorption systems of the quasar in detail. SSB found three C iv systems at z a bs — 2.4509, 2.4742, and 2.8758, 
and one Mg 11 system at z a bs — 0.9598. There is also an LLS at z aos = 2.651, but this system does not have 
associated heavy element lines. Our spectrum ranges from 3800 A to 6190 A. Both Lya and Ly/3 are detected at 
zabs = 2.705 - 2.880. 

z a bs = 2.7719 — There is a weak upward spike at the center of the Ly/3 line in the central trough. Based on this 
feature, we fit the trough with two components having log Nhi = 16.30 and 15.14. We did not detect any metal 
absorption lines. 

z a bs = 2.8113 — Because the central trough has an asymmetrical feature, wc used two components to fit the 
profile. The narrow line at 4625 A was identified as Si 11 A1193 by Chen et al. (1981). However, we identify it as 
Mg 11 A2796 at z abs — 0.654, as the corresponding Mg 11 line of this doublet, Mg 11 A2803, is detected at 4627 A in 
our spectrum. 
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Q0831+1248 (z em = 2.734). — SBS found an absorption system at z abs = 2.0844 with a C iv doublet and Al 11 A1671 
line. Lanzetta et al. (1991) found another system at z ab s = 2.796 with Si 11 and C 11 lines. Our spectrum covers 
from 3790 A to 6190 A. Both Lya and Ly/3 are detected at z abs = 2.695 - 2.734. 

Zabs = 2. 1300 — Metal lines such as C iv and Si iv were detected in this system. As the system is just blueward of 
the emission redshift of the quasar at z em = 2.734, there are few absorption features redward of this system. This 
system may be intrinsically associated with the quasar itself, as the velocity distance from the quasar is small. 

HE0940-1050 (z em — 3.080). — This quasar was discovered during the course of the Hamburg/CfA Bright Quasar 
Survey Reimers et al. (1995) detected four heavy-element systems at z a bs — 2.82 (C iv), 2.32 (C iv), 1.918 (Fe 11, 
Al 11, and Si 11), and 1.06 (Mg 11) with the 4 A resolution spectrum taken with the ESO 1.5m telescope. Reimers et 
al. (1995) did not detect any flux from the quasar below 3200A in a spectrum from IUE, probably due to a Lyman 
limit systems at z a b s <~ 2.82 and/or 2.32. Our spectrum ranges from 3610 A to 6030 A. Both Lya and Ly/3 are 
detected at z abs = 2.519 - 3.080. 

z a bs — 2.8283 — Our spectrum covers Lya to Ly<5 in this system. However, Ly7 and LyS are in low S/N regions. 
Three strong H 1 components with log Nhi = 15.9, 16.4, and 16.0 are located at Av = —300, 0, and +490 km s _1 
from the main component. Seven C iv doublets and five Si iv doublets, C 11 and Si 111 lines were detected in the 
system. 

Zabs = 2.8610 — We fit the central trough with one component having log Nhi = 17.06. However, this component 
may be resolved into multiple narrower components, because the ratio of N\ to N 2 is unusually large, Ni/N 2 ~ 
300. We also detected Si iv and C iv doublets in the system. 

z a bs = 2.9174 — We fit the main trough with five components, referring not only to the profiles of higher orders 
but also to the weak upward spikes in the damped region of the Lya profile. We also detected Si 111 and Si iv lines 
at the redshift of the main component. 

Zabs = 3.0387 — This system is within 5000 km s _1 of the emission redshift of the quasar. At the bluer region 
of the main component there are several unidentified metal lines. They could be identified as C iv lines at z a bs ~ 
2.17, but one of them, at z a b s = 2.158, has a Doppler parameter of b = 28 km s _1 which is unusually large for an 
ordinary C iv line. 

Q1009+2956 (z em =2.644)- — Buries & Tytler (1998b) presented a measurement of the D/H ratio in the metal-poor 
absorption system at z a b s = 2.504. They estimated the D/H ratio in the system to be log(D/H) = — 4.40+Q;gg at 
the 67% confidence level. Our spectrum ranges from 3090 A to 4620 A. Both Lya and Ly/3 are detected at z a b s = 
2.013 - 2.644. 

z a bs = 2.1432 — If the main trough is fit with a single component, it is found to be very broad, b ~ 40 km s . 
We therefore used three components with narrower profiles, 6 = 26, 28, and 31 km s _1 . The fitting model is still 
uncertain, however, as we could refer to only Lya and Ly/3, and both of them are entirely damped. 

Zabs = 2. 4069 — We referred to Lya, Ly/3, Ly7, and hy8 for the profile fitting. The lines of higher orders are 
blanketed by the Lyman continuum of the LLS at z a bs — 2.50. We fit the central trough with one component 
having log Nhi = 18.8, and b = 48 km s _1 . Although this component, with its rather large b value, may be 
resolved into several narrow components, there is too little information to be able to separate them out. This 
system is accompanied by low-ionization metal lines such as O 1, C 11, Si 11, and Si m. 

z a bs = 2.4292 — This system is only 1000 km s _1 to the red of the system at z a bs — 2.407. The spectrum covers 
the system from Lya to Ly6, but they are all damped. We fit the central trough with one component having 
log Nhi — 17.43. However, the ratio of Ni to N 2 is unusually large, ~ 350. The main component may be resolved 
into multiple narrower components. We also detected low-ionization metal lines such as C 11, Si 11, and Si m. 

z a bs — 2.5037 — This system is known to contain a candidate D 1 line. Detailed fitting was applied by Buries & 
Tytler (1998b), and they measured the D/H ratio to be log(D/H) = -4.40t°;og. The H 1 column density of the 
system was estimated by Buries & Tytler (1998b) to be log Nhi = 17.35, which is similar to our result, log Nhi = 
17.26. The D 1 candidate is also detected in our spectrum. Additionally, we detect Si 111 lines. 

Q1017+1055 (z em = 3.156). — SBS found the BAL type features for C iv and Si iv lines at z abs = 2.9720, though 
there were no BAL type Lya absorption lines. Three C iv systems at z a b s — 2.5401, 2.9970, and 3.1101, two Mg 11 
systems at z a bs = 0.974 and 1.2401, and an LLS at z a bs — 3.048 were detected. Our spectrum ranges from 3890 A 
to 6300 A. Both Lya and Ly/3 are detected at z abs = 2.792 - 3.156. 

z a bs — 2.9403 — This system is detected based on only Lya and Ly/3 lines in the regions with a low S/N ratio 
(S/N = 12 at Lya and 8.2 at Ly/3). No metal lines were detected in the system. 

Zabs = 3.0096 — This system is probably identical to the C iv system detected at z a bs — 2.997 in SBS. We detected 
C iv and Si iv doublets in this system, as previously described in SBS. 
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z a bs — 3.054-8 — The central trough was fit with four components having column densities log Nhi — 15.4, 17.1, 
14.9, and 14.0. This model, however, is somewhat uncertain, as the S/N ratio is very low (S/N = 25 at Lya and 
8.9 at Ly/3). The system is accompanied by three weak Si iv doublets. 

z a bs = 3.1120 — This system is detected based on the C iv doublet at z a b s = 3.11 described in SBS, though 
the C iv doublet itself is not covered by our spectrum. We detect two Si iv doublets in the system. Due to the 
asymmetrical feature present in the LyS profile of the central trough, we fit it with two components having log Nhi 
= 15.3 and 15.0. This system is within 5000 km s _1 of the quasar. 

Q1055+4-611 (z a bs = 4-118). — The system at z a b s = 3.32 has been identified as a DLA with logNni = 20.34 (Lu, 
Sargent & Barlow 1998). Storrie-Lombardi & Wolfe (2000) also found another DLA system at z a b s — 3.05, and 
estimated the column density of the DLA to be log Nhi = 20.3; this was confirmed by Peroux et al. (2001). There 
is also an LLS at z aos = 2.90 with large optical depth, t = 4.4 (Peroux et al. 2001). Our spectrum covers a range 
of 4450 A to 6900 A. Both Lya and Ly/3 are detected at z abs = 3.338 - 4.118. 

Zabs — 3.8252 — The central trough was fit with three components with \ogNni = 14.9, 16.0, and 15.6 at Av = 
—65, 0, and 100 km s _1 from the main component. There may exist additional components between Av = and 
+ 100 km s , because there is an excessive residual flux at Lyy and LyS. The system is accompanied by three 
Si iv doublets. The H I lines of orders higher than LyS are all blanketed by the Lyman continuum of the LLS at 

Zabs — 3.93. 

Zabs = 3.8495 — This system is shifted redward of the system at z a bs = 3.82 by just 1000 km s _1 . We fit the 
central trough with two components, using the profiles of Ly<$, Lye, and Ly8 as a reference. The system contains 
Si in and Si iv lines. The Lyman break of the system is not seen, as it is blanketed by the Lyman continuum of 
the LLS at z a bs = 3.93. 

Zabs = 3.9343 — The optical depth of the system at the Lyman break was estimated to be r > 3, which corresponds 
to a column density of log Nhi > 17.7. Using this result as a reference, we fit the main trough with two components 

having a total column density of log Nhi = 17.34. Although there is a data defect at Av = —750 650 km s _1 

from the main component in Lya window, the lines there can be fit using the profiles of higher orders as a reference. 

HS1103+6416 (z em = 2.191). — This quasar was discovered during the Hamburg Quasar Survey (Hagen et al. 1995). 
Kohler et al. (1999) have studied the complex LLS at z a bs — 1-892, using both ultraviolet spectra taken with 
the HST and optical high resolution spectra taken with Keck. They found that the complex absorption lines are 
distributed with a velocity width of ~ 200 km s _1 ; they also found the system to contain at least 11 narrow 
components with various ionization levels. Our spectrum ranges from 3180 A to 5790 A. Both Lya and Ly/3 are 
detected at z a6s = 2.100 - 2.191. 

Q1107+4847 (z em = 3.000). — Carballo et al. (1995) first observed this quasar with a moderate resolution of 40 - 
120 km s -1 , and detected three C iv systems at z aos = 2.697, 2.724, and 2.760. However, two of them, at z aos = 
2.697 and 2.724, are doubtful, because the C iv A1551 components are absent. Our spectrum ranges from 3730 A 
to 6170 A. Both Lya and Ly/3 are detected at z a bs — 2.636 - 3.000. 

z a bs = 2.7243 — The main component has a rather large Dopplcr parameter, b = 48 km s _1 . However, we cannot 
resolve the component into narrower components, as our spectrum covers only Lya and Ly/3 lines, which are both 
strongly damped. This system has four C iv and three Si iv doublets. There are many unidentified metal lines at 
Av = 700 - 1000 km s _1 from the main component in Lya window. Misawa et al. (2007) identified this system as 
a quasar intrinsic system. 

Zabs = 2.7629 — This partial DLA system, with column density log Nhi — 19.13, is accompanied by various metal 
lines such as C iv, Si 11, Si 111, Si iv, and O 1. These metal lines result in complex structure, as described by Carballo 
et al. (1995). Our spectrum covers only Lya and Ly/3 lines, which are both strongly damped. Therefore we used 
the velocity distribution of metal lines as a reference for fitting the H 1 components in the central trough. 

Zabs — 2.8703 — Based on the profiles of Lya, Ly/3, and Ly7 lines, it was possible to separate the components in 
the central trough sufficiently A weak Si iv doublet is also detected in the system. However, we did not detect the 
corresponding C iv doublet, because it is located at the region of the spectral gap. 

Q1157+3143 (z em — 2.992). — There are two LLS candidates at z aos ~ 2.94 and 2.77 in the spectrum of this quasar 
(Kirkman & Tytler 1999). However, the H 1 column density of the LLS at z aos ~ 2.77 has not been evaluated 
exactly, because the LLS at z a bs ~ 2.94 blots out the spectrum below 3600 A, which prevents the detection of H 1 
lines of orders higher than Lyy for the LLS at z aos ~ 2.77. Nonetheless, the system probably has a large column 
density because the system is accompanied not only by high ionization ions such as Si iv, C iv and O vi, but also 
by low ionization ions such as C 11 and Si 11. These low ionization lines are expected to be present in systems with 
large column densities. Our spectrum covers from 3790 A to 6190 A. Both Lya and Ly/3 are detected at z aos = 
2.695 - 2.992. 

Zabs = 2. 7710 — This is one of the two LLSs detected in the direction of this quasar by Kirkman & Tytler (1999). 
As our spectrum covers only Lya and Ly/3 and both are damped, we fit the central trough with one component 
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having log Nhi = 17.63. The system has 7 C iv, 5 Si iv, 2 C n, and 6 Si in lines. The Lyman break is blanketed 
by the Lyman continuum region of another LLS at z abs ~ 2.94. 

z a bs = 2.8757 — Around the center of this system there are two H I lines with very similar column densities, log Nhi 
= 15.54 and 15.66, and with a velocity separation of 210 km s _1 . Interestingly, only one of them is accompanied 
by various metal lines, such as Si in, Si iv, and C iv. 

z abs = 2.9437 — This system is another LLS detected by Kirkman & Tytler (1999) with log N H i = 17.44. This 
LLS has a companion H I line with similar column density of log Nhi = 17.16 having a separation of 250 km s _1 . 
Si iv, and Si 11 lines were detected in the system. Two O 1 lines were also detected. However, the detection of O 1 
lines is tentative, because they have shifted as much as 400 km s _1 blueward of the main component. 

Q1208+1011 (z em = 3.803). — This is a candidate gravitationally lensed quasar (Bahcall ct al. 1992). Two point 
sources with magnitudes of V = 18.3 and V = 19.8 are separated by 0".47. The spectra of both objects have strong 
Lya+N v and O vi+Ly/3 emission lines at z em ~ 3.8, and many common absorption features. The low resolution 
(25 A) spectrum of Bahcall et al. (1992), however, could not resolve whether these objects are (i) gravitational 
lensed images, or (ii) two quasars with a projected separation of a few kpc. Our spectrum ranges from 3730 A to 
6170 A. Both Lya and Ly/3 were detected at z abs = 2.636 - 3.803. 

z a bs — 3.3846 — The spectrum covers all Lyman series orders of this system, but they are positioned in regions 
with low S/N ratio (e.g., S/N = 24 at Lya). As the central trough is not separated into narrow components up to 
LylO, we fit the trough with a single component having log Nhi = 17.35. We did not detect any metal lines in the 
system. 

Zabs = 3.4596 — The Lya lines of this system are blended with C iv lines at z a bs ~ 2.50. Our fitting model, when 
applied to the Lya, Ly7, and Ly5 lines, does not agree with the wings of Ly/3, which suggests that the wings are 
not Ly/3 lines at z a bs ~ 3.46 but other, lower-redshift lines such as Lya. 

Zabs = 3.5195 — We fit the central trough with two components based on the features of Ly8 and LylO. The 
profiles of Lya - Ly7 are all saturated. 

Zabs = 3.7206 — We fit the system with a component with a rather small column density, log Nhi = 15.48. No 
metal lines associated with the system were detected in our spectrum. 

Q1244+H29 (z em = 2.960). — This quasar was discovered in the Hamburg/CfA Bright Quasar Survey. No results of 
detailed analysis for this quasar have been published. Our spectrum ranges from 3370 A to 4880 A. Both Lya and 
Ly/3 were detected at z a bs — 2.285 - 2.960. 

Zabs — 2.9318 — The central trough, which has an asymmetrical profile, was fit with four components using higher 
order H I lines as a reference. No metal lines were detected in the system. The system is within 5000 km s _1 of 
the quasar. 

Q 1251 +3644 ( z em = 2.988). — There is a Lyman break at 3303 A in this system, produced by the LLS at z abs = 2.614 
(Stenglcr-Larrea et al. 1995). Our spectrum ranges from 3790 A to 6190 A. Both Lya and Ly/3 are detected at 
zabs = 2.695 - 2.988. 

Zabs = 2.8684 — Only Lya and Ly/3 lines are covered by the spectrum. Although both orders are affected by 
several data defects, the fit model is good. We did not detect any metal lines in this system. 

Q1330+0108 (z em = 3.510). — A Lyman break at 4034 A, produced by the LLS absorber at z a6s = 3.414, was detected 
(Stcngler-Larrca ct al. 1995). Our spectrum ranges from 4030 A to 6450 A. Both Lya and Ly/3 are detected at 
zabs = 2.929 - 3.510. 

Q1334-0033 (z em = 2.801). — Although Wolfe et al. (1995) attempted to find DLA systems in the spectrum of this 
quasar, they found no candidates for DLA systems in the observed redshift range, 1.634 < z < 2.745. No absorption 
systems in the spectrum of this LBQS quasar have been published. Our spectrum ranges from 3730 A to 6170 A. 
Both Lya and Ly/3 are detected at z a bs — 2.636 - 2.801. 

z a bs — 2.7572 — Our spectrum covers only Lya and Ly/3 lines, and the S/N ratio around Ly/3 is very low, S/N < 
10. Therefore the fitting model is not good. However, the existence of three C iv doublets in the system strongly 
suggests that the system is real. This system is within 5000 km s _1 of the quasar. 

Q1337+2832 (z em = 2.537). — This quasar was discovered by the CFHT/MMT survey (Crampton ct al. 1988,1989). 
Our spectrum covers a range of 3170 A to 4710 A. Both Lya and Ly/3 are detected at z a b s — 2.091 - 2.537. 

Zabs = 2.4336 — The optical depth of the system at the Lyman break was estimated to be r > 3, corresponding 
to a column density of log Nhi > 17.7. In our spectrum, we used only lines of orders up to LyS, because the S/N 
ratios of the higher orders are very low (S/N = 2.5 at Lye). We fit the central trough with four components having 
a total column density of log Nhi = 18.9, based on the profile of Ly<5. While high ionization lines such as the C iv 
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and Si iv doublets of this system are not covered by the spectrum, low ionization lines (e.g., C n, Si n, Si in, and 
O i) were detected. 

z a bs — 2.5228 — This system is within 5000 km s _1 of the quasar. The central trough was fit with two components 
having column densities of log Nhi = 15.81 and 14.44, because the trough has an asymmetrical profile at Lya and 
Ly/3. We did not detect any metal lines in the system. 

QU22+2309 (z em = 3.611). — Kirkman & Tytler (1997b) have detected the LLS at z abs = 3.3816 with high ionization 
metal lines such as C iv and O vi. The H I column density of the LLS was estimated to be logNni — 20.6 assuming 
the LLS is collisionally ionized, and log Nhi — 19.9 if it is photoionized. Our spectrum ranges from 3740 A to 
6180 A. Both Lya and Ly/3 arc detected at z abs = 2.646 - 3.611. 

z a bs — 3.3825 — The spectrum covers the wavelength range from Lya to the Lyman limit of this system. The 
S/N ratio is high over all orders. We fit the central trough with two components having similar column densities, 
log Nhi = 16.5 and 16.3, using the profile of orders higher than Lye as a reference. 

z a bs — 3.5362 — The wavelength range of our spectrum also covers all Lyman orders. As the Lya lines of the 
system are located in the gap of the spectrum, we fit them by reference to the higher orders. We resolved the 
central trough into 6 narrow components with column densities of log Nhi = 14.5 - 15.9. Three C n and four Si in 
lines were also detected. This system is within 5000 km s _1 of the quasar. 

Q1425+6039 (z em — 3.165). — This quasar is one of four very luminous quasars at z em > 2 found by the Second Quasar 
Survey (Stepanian et al. 1990 and references therein). Stepanian et al. (1996) carried out follow-up spectroscopy 
for this quasar with the 6m telescope of the Special Astrophysical Observatory, and detected a strong DLA system 
with log Nhi = 20.4 at z abs = 2.826. This system contains not only low-ionization ions (e.g., C 11, Si 11, Fe 11, and 
Al 11), but also high- ionization ions (e.g., O vi, N v, Si iv, and C iv). Our spectrum ranges from 3730 A to 6170 A. 
Both Lya and Ly/3 are detected at z a bs = 2.636 - 3.165. 

Zabs = 2.1100 — This partial DLA system has a large column density of log Nhi = 19.4, and results in strong 
damping wings, preventing us from detecting the weak components around the DLA system. We fit the central 
trough with seven components. Three H 1 components near the center of the system are located at the same redshift 
as three Si 111 components in the system. We fit the Doppler wings with four components using the profiles of Ly/3 
lines as a reference. There are also Si 11, Si iv, and C iv lines in the system. Although this system was classified as 
a possible candidate for a quasar intrinsic system, it was less reliable because of this DLA-like structure (Misawa 
et al. 2007). 

Zabs = 2.8258 — This DLA system has a large column density of log Nhi = 20.0, which results in strong damping 
wings. Almost all of the H 1 lines within 1000 km s _1 of the main component are blanketed by the wings. This 
system is not suitable for our study of the structure of H 1 systems. In this DLA system, 7 C iv, 7 C 11, 9 Si iv, 2 
Si in and 6 Si 11 lines were detected. 

Zabs — 3.0611 — We could not detect the Lyman break of the system around 3740 A due to a low S/N ratio. 
The Lyi5 profile has two upward spikes on the bottom of the central trough. If the trough is fit with three narrow 
components based on these spikes, the Doppler parameters of them are found to be b ~ 7 km s _1 , which is too low 
for H 1 lines. Therefore we fit the trough with two components, using the asymmetrical Lya profile as a reference. 

Zabs = 3.1356 — The central trough was fit with two components having log Nhi — 16.7 and 16.2, using the 
profiles of orders higher than Lye as a reference. The region between Av = —500 and —300 km s _1 from the main 
component in Lya window is affected by the gap of the spectrum. Therefore we fit the region based on lines of 
order higher than Ly/3. The system is within 5000 km s _1 of the emission redshift of the quasar. 

QU42+2931 (z em = 2.610). — Carballo et al. (1995) found three C iv systems at z abs = 2.330, 2.439, and 2.474 in 
the spectrum of this quasar. The strongest system at z abs = 2.439 was found to be accompanied by low-ionization 
lines such as O vi and C 11. Carballo et al. (1995) also detected two Lya absorption lines with log Nhi > 16 at z abs 
= 2.555 and 2.617, though no metal lines were detected in their redshift range. Our spectrum ranges from 3740 A 
to 6180 A. Both Lya and Ly/3 are detected at z abs = 2.646 - 2.670. 

Q1526+6101 (z em = 3.020). — This quasar was discovered during observations taken with the NRAO Green Bank 300 
foot telescope (Becker et al. 1991). Storrie-Lombardi & Wolfe (2000) confirmed that there were no candidate DLA 
systems in their observed redshift range, 1.955 < z < 2.980. Our spectrum ranges from 3460 A to 4980 A. Both 
Lya and Ly/3 are detected at z abs — 2.373 - 3.020. 

Zabs = 2.9151 — We detected the system in a region of low S/N ratio (e.g., S/N = 24 at Lya, and 7.5 at Ly/3). We 
fit the central trough with three components using the profile of Ly/3 as a reference. No metal lines were detected 
in the system. The system is within 5000 km s _1 of the quasar. 

Q1548+0911 (z em = 2.149). — SBS found a fairly weak, but unambiguous, Mg 11 doublet at z abs = 0.7708 in the 
spectrum of the quasar. They also detected three C iv doublets at z abs = 2.2484, 2.3195, and 2.4915. Our 
spectrum ranges from 3730 A to 6180 A, which does not cover the heavy element systems detected in SBS. 
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Q1554+3749 (z em = 2.664). — This quasar was discovered by the Palomar Transit Grism Survey (PTGS; Schneider, 
Schmidt, & Gunn 1994). No absorption systems in the quasar have been published. Our spectrum ranges from 
3240 A to 4770 A. Both Lya and Ly/3 are detected at z abs = 2.159 - 2.664. 

z a bs = 2.6127 — The central trough was fit with two components having log TV hi = 18.0 and 14.5. Although the 
trough may be separated into more than two components, we cannot resolve it with our low S/N spectrum. The 
strong lines at Av = —450 km s _1 and —250 km s _1 from Lya of the main component are not Lya lines, because 
there are no corresponding Ly/3 lines. They are also not Ly/3 lines at higher redshift, since the corresponding 
redshift, z a b s ~ 3.27, is higher than the emission redshift of the quasar, z em = 2.664. Therefore they are likely to 
be very strong metal absorption lines. 

HS1700+6416 (z em = 2.722). — This quasar has been well studied by both ground and space telescopes (Reimers et al. 
1989; Sanz et al. 1993; Reimers et al. 1992; Rodrfguez-Pascual et al. 1995), as 6 (sub-)LLS candidates have been 
identified at z abs = 2.4336, 2.1681, 1.8465, 1.725, 1.1572 and 0.8642 with 16.0 < logN HI < 18.3. There are also 
4 C iv doublets at z a b s = 2.7444, 2.7102, 2.5784 and 2.308. One of them is at z a \, s > z em . Our spectrum ranges 
from 3730 A to 6180 A. Both Lya and Ly/3 are detected at z a b s = 2.636 - 2.722. We detected the C iv systems at 
z a bs = 2.7 '444 and 2.7102 in our spectrum. However, the corresponding H 1 lines are very weak, log Nhi ~ 14.0. 

Q1759+7529 (z em = 3.050). — Outram et al. (1999) studied this quasar in detail. They detected a DLA system at 
z a bs — 2.625 with log Nhi = 20.76. We also found a partial DLA system at z a b s — 2.910 with log Nhi = 19.80, 
which is accompanied by an LLS with logNni = 17.02 having a separation of 420 km s . In the spectrum of the 
quasar, 9 C iv systems were detected at z abs = 1.8848, 1.935, 2.4390, 2.484, 2.7871, 2.795, 2.835, 2.84, and 2.896. 
Three systems at z a bs = 1-935, 2.484 and 2.625 have complex structures with both low- and high-ionization lines. 
Galactic absorption lines of Na 1 AA5892,5898 were also observed. Our spectrum ranges from 3580 A to 5050 A. 
Both Lya and Ly/3 are detected at z a bs — 2.490 - 3.050. 

Zabs = 2.7953 — We fit the central trough with four components having column densities of log Nhi = 13.4, 14.0, 
15.3, and 14.9, based on the profiles of Ly7 and Ly<5. There are no metal lines in the system. 

Zabs = 2.8493 — The central trough was fit with five components referring to the profile of LyS. Both sides of 
the Ly/3 line could not be fit well by our model. If additional components were included in the fit, the model was 
found to absorb too much at Lya as compared with the observed spectrum. No metal lines were detected in our 
spectrum, though the corresponding C iv and Si iv lines are not in the observed range. 

Zabs = 2.9105 — The observed spectrum could not be correctly normalized, as this partial DLA system has strong 
damping wings. We fit the strong trough with 8 components, using the profiles of lines of higher order than Ly/3. 
The column density of the main component is logNni = 19.90, which is almost identical to the value determined 
by Outram et al. (1999), \ogN H i = 19.80. 

Q1937-1009 (z em = 3.806). — LLS with \ogN HI = 17.86 was detected at z abs = 3.572 by Buries & Tytler (1997). The 
LLS is an ideal system for obtaining an estimate of the primordial value of D/H, because the system has very low 
metallicity - less than 10 -2 solar abundance. The D/H ratio has been evaluated to be D/H = 3.3±0.3 x 10~ 5 at 
the 67 % confidence level (Tytler, Fan, & Buries 1996; Buries & Tytler 1998a). The LLS is accompanied by various 
metal lines, such as C 11, C 111, C iv, N 111, Si 11, Si 111, Si iv, Fe 11, and Fe m. Our spectrum ranges from 3890 A to 
7450 A. Both Lya and Ly/3 are detected at z abs = 2.792 - 3.806. 

Zabs — 3.5725 — Buries & Tytler (1998a) found a candidate D 1 line in this system. They measured the column 
density of the H 1 line to be logNni — 17.86. We fit the central trough with two components having column 
densities of log Nhi = 17.94 and 15.89. Various metal lines such as C 11, C 111, C iv, Si 11, Si in, and Si iv were also 
detected in the system, though all of them have only a single component. 

HS1946+ 7658 (z ern — 3.051). — This quasar was discovered by Hagen et al. (1992) using objective prism observations 
at the 80 cm Calar Alto Schmidt telescope. Three metal absorption systems have been found. Two of them are 
highly ionized systems at z a bs — 3.049 and 2.843 with C iv, Si iv, and N v lines. The other one is a low ionization 
system at z a b s — 1.738 with Mg 11 and Fe 11 lines (Hagen et al. 1992; Sadakane et al. 1993). Our spectrum ranges 
from 3890 A to 6300 A. Both Lya and Ly/3 are detected at z a bs = 2.792 - 3.051. The system at z a b s = 2.843 is 
located in the spectral gap. 

z a bs = 3.0498 — As this system is within 1000 km s _1 of the emission redshift of the quasar, there are only a few 
absorption lines redward of the main component. The detection of N v lines in this system strongly suggests that 
the system is highly ionized by the UV flux of the quasar. This system is probably intrinsically associated with the 
quasar itself, because Misawa et al. (2007) found partial coverage in the N v doublet. Other metal lines, such as 
C iv, Si 11, Si in, and Si iv, were also detected. 

Q2223+2024 (z em = 3.560). — This quasar was first discovered by the MIT-Green Bank III 5 GHz survey (Griffith et 
al. 1990). Storrie-Lombardi & Wolfe (2000) confirmed that there was no DLA candidate in their spectrum with a 
range of 2.101 < z < 3.514. No detections of absorption systems have been published for this quasar. Our spectrum 
has a range of 4120 A to 6520 A. Both Lya and Ly/3 are detected at z a bs = 3.017 - 3.560. 
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Q2344 +2024 (z em = 2.763). — This quasar is only ~ 5' separated on the sky from the quasar Q2343+125. Q2344+2024 
has C iv systems at z a b s = 2.427 and 2.429, while Q2343+125 has corresponding C iv systems at z a b s — 2.429 
and 2.431 (SBS). The velocity separation of these systems along the line of sight is less than 1000 km s _1 . These 
common C iv systems could be produced by a Mpc-scale absorber, such as a cluster of galaxies at the given redshift, 
although Misawa et al. (2006) did not find any Ha emitting galaxies at z ~ 2.43 in this pair quasar field, down to 
f(Ra) of 1.6xl0" 17 erg s" 1 cm" 2 . Q2344+2024 has other C iv systems at z abs = 2.2754, 2.4265, 2.4371, 2.6964, 
2.7017, and 2.7817. Our spectrum has a range of 3410 A to 4940 A. Both Lya and Ly/3 were detected at z a bs = 
2.324 - 2.763. 

Zabs = 2. 4261 — The spectrum covers only Lya and Ly/3 in this system. The model fit is poor, because both Lya 
and Ly/3 are in regions with low S/N ratio (S/N = 16 at Lya and 6.8 at Ly/3). Weak Si iv doublet lines are also 
detected in this system. 

Zabs = 2.6356 — We detected only a few H I components in the system; this is partially due to low S/N ratio. 
The Doppler parameter of the main component, b = 62 km s _1 , is larger than would be expected for a H I line, 
suggesting that the main component may be separated into several narrow components. We did not detect any 
metal lines in the system. 

z a bs = 2.7107 — The component that has the second largest column density, log Njji = 15.26, is separated by 
—750 km s^ 1 from the main component, which has log Njji — 16.64. It is unclear whether the two components are 
physically related to each other. This system is within 5000 km s _1 of the emission redshift of the quasar. 

z a bs = 2.7469 — The H I component that has second largest column density, log Njji — 16.27, is rcdshiftcd 
+560 km s _1 from the main component, which has log Njji — 16.67. It is unclear whether the two components 
are physically related. There are no corresponding metal lines in the system. This system is within 1000 km s _1 
of the emission redshift of the quasar. 
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Table A10 

(Sample) Absorption lines in 86 H i systems 



(1) (2) (3) (4) (5) (6) (7) (8) (9) 

No. \ obs Av z abs log iV o-(logJV) b cr(b) ID 

(A) (kms- 1 ) (cm- 2 ) (cm- 2 ) (kms" 1 ) (kms" 1 ) 



Q0004+1711 (z,- m =2.890j 
z at , s =2.8284 



1 


4639.9 


-917.4 


2.81671 


13.505 


0.018 


29.48 


1.43 


H I A1216 


2 


4641.1 


-837.6 


2.81772 


13.246 


0.027 


27.24 


1.99 


H I A1216 


3 


4643.0 


-717.1 


2.81926 


12.739 


0.069 


11.86 


1.81 


M I 


4 


4644.0 


-650.4 


2.82011 


13.048 


0.164 


41.21 


19.08 


H I A1216 


5 


4645.8 


-531.6 


2.82162 


14.459 


0.019 


28.88 


0.47 


H I A1216 


6 


4649.6 


-289.2 


2.82471 


13.244 


0.040 


31.27 


3.82 


H I A1216 


7 


4654.1 


0.0 


2.82840 


15.506 


0.040 


43.76 


0.61 


H I A1216 


8 


4656.1 


132.7 


2.83009 


13.087 


0.057 


25.87 


3.96 


H I A1216 


9 


4656.6 


161.0 


2.83045 


12.996 


0.042 


5.13 


0.52 


M I 


10 


4656.8 


177.5 


2.83066 


13.103 


0.035 


8.31 


0.88 


M I 


11 


4657.2 


203.4 


2.83100 


13.136 


0.046 


31.44 


3.96 


H I A1216 


12 


4667.6 


869.2 


2.83951 


13.722 


0.014 


73.98 


3.36 


H I A1216 


13 


4669.6 


999.2 


2.84117 


14.130 


0.030 


22.09 


0.78 


H I A1216 



Table All 

(Sample) Detected Metal Absorption lines in 86 H i systems 



(1) (2) (3) (4) (5) (6) (7) (8) (9) 

quasar z sys Av z abs log N c(log N) b ^ity ID 

(km s" 1 ) (cm -2 ) (cm" 2 ) (km s" 1 ) (km s" 1 ) 



Q0001 + 1711 


2.8707 


-26.0 


2.87034 


14.052 


0.083 


26.77 


2.00 


C in 






-23.9 


2.87037 


14.673 


0.038 


21.05 


0.82 


C II 






-23.7 


2.87037 


14.063 


0.037 


17.27 


0.39 


Si II 


Q0014+8118 


2.7989 


-86.1 


2.79782 


12.634 


0.015 


21.72 


0.93 


Si in 






-57.7 


2.79818 


12.387 


0.021 


6.60 


0.26 


Si in 






-56.5 


2.79820 


12.950 


0.016 


6.96 


0.27 


C iv 






-56.2 


2.79820 


12.537 


0.010 


6.91 


0.17 


Si iv 






14.9 


2.79910 


13.032 


0.024 


18.89 


1.14 


C iv 






17.5 


2.79913 


12.204 


0.082 


11.78 


2.55 


Si iv 






34.4 


2.79935 


12.055 


0.206 


5.72 


2.02 


C iv 






34.7 


2.79935 


11.456 


0.312 


5.19 


2.23 


Si iv 






58.2 


2.79965 


12.536 


0.046 


14.24 


1.05 


Si iv 






61.7 


2.79969 


13.525 


0.015 


32.73 


1.20 


Si in 






67.3 


2.79977 


13.576 


0.022 


25.91 


1.59 


C iv 
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Fig. All. — (Sample) Velocity maps of the lowest five orders of Lyman series (i.e., Lya, Ly/3, Ly7, Ly5, and Lye) in a H I system at z a t, 3 
= 2.8284 in the spectrum of Q0004+1711. A histogram just above the zero flux is la error spectrum. Tick marks above the flux spectrum 
denote the positions of H I and M I lines. Metal lines in the velocity window are not marked because they are not physically related to the 
H I system. 



